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THE REACTION! 


ABSTRACT 


The reaction has been observed using thin targets tritium absorbed 
zirconium bombarded singly charged helium ions energy from 0.5 Mev 
1.9 Mev. The cross section rises smoothly with energy fashion charac- 
teristic direct radiative capture The ratio the intensities 
transitions the first excited state 478 kev and the ground state 0.40+ 
0.05. The angular distribution is, within errors, isotropic 0.8 Mev 
but significantly higher than 90° for 1.6 Mev. The total capture 
cross section mean alpha-particle energy the target 1.32 Mev 
3.58+0.60 microbarns and the corresponding astrophysical factor 0.064+ 
0.016 kev barn the center mass system. The results are compared with 
recent theoretical results the direct radiative capture process. 


INTRODUCTION 


the last few years there has been considerable interest direct reactions 
which take place without the formation compound nucleus and which, 
many cases, depend the nuclear wave function and even beyond the 
conventional nuclear radius. particular class these reactions the direct 
radiative capture process. This process contributes significant gamma-ray 
yield the energy regions between resonances for number reactions 
low bombarding energy, such (Warren al. 1954), 
(Tanner 1959), and T(a,y)Li’ (Riley 1958 and Griffiths 1959). 

have extended the earlier measurements made this laboratory 
the reaction. The reaction has also been studied Holmgren and 
Johnston (1959), who studied the mirror reaction well. The 
cross sections are over 2000 times greater than earlier estimates given 
Cameron (1957) based work Salpeter (1952), which did not consider 
the direct reaction mechanism. For this reason these reactions, particularly 
the mirror reaction, are importance nucleogenesis since they provide 
second method bypassing the mass nuclei pointed out Fowler (1958) 
and Cameron (1958). addition, the results are interest because recent 
developments direct capture theory Christy and Duck (1961) and 
Tombrello and Phillips (1961). The latter authors have interpreted the 
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radiative capture terms “two-body model for and Be’ 
which assumes that low excitation the mass nuclei can thought 


EXPERIMENTAL 


Targets tritium absorbed thin evaporated layers zirconium 
platinum backings, kindly supplied the Oak Ridge National Laboratories, 
were bombarded with microamperes singly charged helium ions from 
the University British Columbia Van Graaff generator the energy 
range from 0.5 2.0 Mev. The capture gamma rays were observed with well- 
shielded sodium iodide scintillation counters. Because the low gamma-ray 
yield, backgrounds, particularly that from were reduced 
thorough cleaning parts struck the beam and placing liquid air 
trap the beam tube just front the target holder. 

Scintillation Counter Efficiency 

The counter used most this work consisted in. long in. 
diameter sodium iodide crystal mounted Dumont K1213 photomultiplier 
operated 900 volts, with pulse height analysis provided 100 channel 
Computing Devices Canada kicksorter. Five inches 
were used the counter reduce the background from the generator and 
from radiothorium and the concrete walls the laboratory. Provision 
was made rotate the counter assembly through +135° the incident beam 
direction. 

theoretical and experimental study the efficiency this counter for 
y-rays the energy range from 0.5 Mev Mev has been made with the 
results shown Fig. Although the gamma rays lie the energy 


Efficiency 


0 3 4 6 8 10 12 14 16 18 20 22 
Energy- Mev 


Fic. efficiency for the in. diameter in. long Nal counter showing 
the efficiency for bias equal half the y-ray energy, and for bias set just below the full 
energy peak, The dotted curves show the total theoretical primary absorption (upper 
curve) and the part the total theoretical primary absorption which initially leaves more 
than half the y-ray energy the crystal. 
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region from 0.4 3.5 Mev the efficiency Mev was also required because 
the tritium content the targets was measured terms the yield 20-Mev 
gamma rays from the reaction. The efficiency, used here 
defined for given bias the counting system the following relation for 
the number counts above the bias, 


where the first factor the number gamma rays per unit solid angle, 
the area the end face the counter, the distance from the end face 
the gamma-ray source, and the center’’ distance. This distance 
measured from the end the crystal point inside the crystal such that 
the counting rate varies the inverse square the distance the gamma-ray 
source from this point. function gamma-ray energy and obtained 
from the intercept the distance axis plots the square root the 
counting rate versus distance. These plots are remarkably straight lines for 
distances several meters less than the length the crystal. Figure 
shows the efficiency, for bias equal half the y-ray energy used 
most this work, well the efficiency for bias set just below the full 
energy peak, The two lower energy points were obtained terms 
absolutely calibrated sources and The 6.14-Mev point was 
determined simultaneous counting alpha particles from the 
reaction with thin window proportional counter suggested Van Allen 
and Smith (1941). The upper the curve was obtained inter- 
comparison with standard brass geiger counter similar that used and 
calibrated Barnes (1952). Details this work will published 
later date. 
Tritium Targets 

The characteristics the four targets used are given Table These 


TABLE 
Target tritium content 


Tritium content, atoms/cm? (kev) 


Ratio, 0.8-Mev 1.5-Mev 
Number B-rays T/Zr protons alphas 
1.73 1.75 1.72 1.0 32.5 190 
1.35 1.30 2.41 1.8 25.5 148 
2.97 3.81 385 


3.53 


consisted weighed thickness zirconium evaporated onto 0.01 in. thick 
platinum disks and saturated with tritium. They were fixed thin copper 
backing mounted cooled copper rod. There was appreciable loss 
yield after much hours running with microamperes 1.5-Mev 
alphas spot. The uniformity tritium surface density was checked 
observing the yield various places the surface and found 
better than 5%. 
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The tritium content was estimated the suppliers terms the zirconium 
thickness assuming one one atomic ratio Zr, and measuring 
the ionization produced the 18-kev beta activity from the target. The 
suppliers indicated that the atomic ratio could vary more than 50% and 
for this reason have made independent estimates the tritium content 
shown under the heading ‘‘by y-rays” Table First, the relative tritium 
contents were obtained measuring the neutron yield from the 
reaction 1.5 Mev using y-ray insensitive plastic phosphor neutron counter 
(Griffiths al. 1959). Second, the tritium was estimated absolutely terms 
the 20-Mev yield from the reaction 800 kev, this energy 
being chosen well below the 873-kev resonance fluorine because the fluorine 
contamination the targets was appreciable. The absolute tritium content 
was based the cross section given Perry and Bame (1955) 
for the mean proton energy the target cm?/steradian 90° 
for target number 45° the beam and mean proton energy 755 kev). 
The final two columns Table give the energy loss for protons and alpha 
particles passing perpendicularly through the target, assuming one one 
atomic ratio and using figures interpolated from the stopping cross-section 
data given Whaling (1958). 

typical 20-Mev y-ray spectrum shown Fig. Since the 


Fluorine 

| 

J 


Half Energy Bios 


ENERGY- MEV 


Fic. Typical spectrum used for estimating the tritium content the targets 
showing the half energy bias and extrapolated tail. 


cross section obtained Perry and Bame (1955) using gas target was based 
similar scintillation counter spectrum obtained with crystal com- 
parable size have extrapolated the tail the spectrum zero pulse height 
shown and assumed that the number counts under this spectrum corre- 
sponded the total theoretical primary absorption the crystal, was done 
Perry and Bame. For these measurements the counter was centimeters 
from the target, therefore only small solid angle corrections were required. 
The gamma-ray flux estimated this method 25% higher than that 
estimated terms the efficiency for half energy bias shown Fig. for 
20-Mev y-rays obtained intercomparison with the standard geiger counter. 
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This suggests that the tail the spectrum shown contains considerably 
more pulses than can accounted for primary absorption events the 
crystal. This discrepancy, although may introduce some uncertainty the 
cross section, should cancelled deducing the tritium content 
the target since have made similar interpretation the scintillation 
counter spectrum for comparable sized counter was done establishing 
the cross section. The ratio tritium content estimated y-rays that 
estimated from the zirconium thickness, assuming one one atomic ratio, 
shown Table The agreement very good for targets and which 
were used for the absolute yield measurements. Target was used for angular 
distribution and energy measurements and target thick target, was used 
for investigating secondary neutron effects discussed below. decay correction 
nearly 10% was made for the tritium decay between the time target was 
made and when was used; not known the He* daughter the decay 
escapes from the target. 


RESULTS 
Gamma-ray Spectra 


Figure shows typical y-ray spectra observed for alpha-particle bombarding 
energies 1.6 Mev and 0.8 Mev with beam currents microamperes 


1000 


a 
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Fic. y-ray spectra two bombarding energies. 


target number The inset shows the energy region interest. the 
energy region from Mev background approximately counts per 
channel has been subtracted from the 1.6-Mev curve, which was obtained 
minutes with integrated current 5000 microcoulombs. this 
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1-Mev energy region the time-dependent background was counts per minute 
and the beam-dependent background was counts per minute 1.6 
Mev and counts per minute 0.8 Mev measured clean piece 
zirconium. The mean alpha-particle energies the target for these runs 
were 1.32 Mevand 0.56 Mev with the target 45° the incident 
beam. 

addition the three y-rays apparent that there sub- 
stantial background from the target particularly low energies. Although the 
presence two overlapping y-ray spectra 3.0 and 2.5 Mev may not 
clear Fig. because the separation energy the two y-rays, and 
almost the same the separation the multiple peaks produced each 
y-ray, the presence the 2.5 Mev, the first excited state Li’ con- 
firmed the presence the 0.477-Mev ground state transition, from the 
first excited state. addition, the spectra not have the shape characteristic 
single y-ray this energy region. order separate the contributions 
the spectra from and standard spectra were obtained from radio- 
thorium (2.62 Mev), (3.09 Mev), Mev), and 
(6.14 Mev). From these, standard spectra could constructed 
for intermediate energies. The peaks the standard spectra had 
broadened slightly from the shapes obtained for monochromatic y-rays 
order that the decomposition procedure would result reasonably smooth 
background, because the spectra are broadened due the appreci- 
able target thickness. The procedure was fit standard spectrum the 
upper photopeak for with allowance for the background above this peak 
and then subtract it. Finally spectrum was subtracted, leaving 
smooth background. decomposition this type shown Fig. for data 
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Fic. y-ray spectrum 1.6 Mev showing the decomposition into 
components and the remaining background (broken line). 
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taken bombarding energy 1.6 Mev showing the residual background 
which comes from the target, since ordinary time-dependent and beam- 
dependent backgrounds have been subtracted. 

The y-ray energy was found vary where the mean alpha 
energy the target. The mean alpha energies calculated from the observed 
y-ray energies agreed better than with those calculated from the target 
thickness data. rough estimate the intensity agreed with the intensity 
obtained after separating the spectra, confirming that these two y-rays 
are cascade. 

Measurements the y-ray energy function angle were made 
0°, 45°, 90°, and 135° for bombarding energy 1.6 Mev. These clearly 
showed the expected Doppler shift amounting 42+5 kev for between 
90° and good agreement with calculated shift kev. 


Branching Ratio 

From the decomposed spectra and knowledge the energy dependence 
the counter efficiency was possible determine the ratio intensities 
the transitions the first excited state and the ground state. Within 
the experimental accuracy the ratio was independent angle 
indicating that and have practically the same angular distribution. The 
branching ratios are 


0.45+0.03 1.6 Mev 1.32 Mev), 


Angular Distributions 
the assumption direct radiative process one would expect the 
reaction low energies proceed electric dipole radiation 


from continuum state wave a-particles) the ground state 


and first excited state. this case the angular distributions should 
isotropic. The preliminary results Riley al. (1958) indicated that 
1.6-Mev bombarding energy the yield was greater than 90°. For this 
reason more detailed angular distribution results have been obtained 1.6 
Mev and 0.8 Mev. Only the ground state radiation was analyzed since the 
branching ratio independent angle within the accuracy about 5%. 

800 kev 0.56 Mev), measurements were made and 90° only, 
giving 90° yield ratio 1.05+0.06. Thus this energy the angular 
distribution isotropic within the errors the measurements. 

1.6 Mev, spectra were obtained for and 45°, 90° and 135° 
each side the target. The target plane was positioned either 90° 
+45° the beam direction, depending the angle observation, and 
corrections were made for target thickness and absorption the backing 
due the change target angle. The results for six different runs each 
angle were averaged give the intensity ratios shown Table II. The errors 
given are based counting statistics only but they are close the root 
mean square errors the individual runs from the mean for each angle. The 
results were symmetrical either side better than 4%. Finally small 
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TABLE 
Angular distribution ratios for 1.6 Mev 


Intensity ratio Corrected 
0°/90° 1.30+0.04 1.28+0.04 
45°/90° 1.08+0.06 1.06+0.06 
135°/90° 1.02+0.06 


corrections were applied for the finite solid angle the detector and for 
center mass motion give the corrected ratios shown Table II. should 
pointed out that the target-dependent background left over after separating 
the standard spectra for and was not itself isotropic; this discussed 
the next section. However, felt that the decomposition procedure used 
for the spectra has effectively removed any effect the background anisotropy 
from the angular distribution for 
Neutron Background 

The target-dependent background, which was subtracted from the data 
taken for the 1.6-Mev angular distributions, was very markedly anisotropic 
with distribution that could fitted the function (1+0.35 cos 
addition, the spectrum this background which remained after separating 
out the standard y-ray spectra, was much like the background produced 
the sodium iodide counter neutrons. was unlikely that the 
reaction contributed appreciably independent checks the 
and reactions indicated very low carbon contamination 
the target. Using thick target (number bombarded 1.6-Mev 
a-particles, the neutron spectrum, yield, and angular distribution were investi- 
gated with y-ray insensitive neutron counter (Griffiths al. 1959). This 
counter has poor energy resolution, but comparison with monoenergetic 
calibration spectra and neutron spectra was concluded that the 
neutrons from the tritium target had energies between and Mev. This 
seems rule out the possibility that the neutrons were produced small 
amount molecular deuterium the mass beam which would have pro- 
duced 14-Mev neutrons the reaction. seems most likely that 
the neutrons observed are result secondary process; tritons scattered 
incident a-particles collide with further target tritium atoms produce 
neutrons similar the secondary neutron production which has 
been observed heavy ice targets under proton bombardment (Singh al. 
1959). The anisotropy the neutron yield from target number 1.6 Mev 
was measured and found agree within 10% with the angular distribution 
the target-dependent background observed the scintillation counter. 
addition, rough check showed that the ratio the neutron production 
1.6 Mev that 0.8 Mev was also agreement with the ratio the 
residual scintillation counter backgrounds these energies. These results then 
support the view that the greatest part the target-dependent background 
due neutrons from the target, likely produced scattered tritons re- 
acting with other target tritons. This secondary neutron effect should even 
greater for gas target where there are additional zirconium atoms 
contribute stopping the recoil tritons. 
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Cross Section 

The absolute 90° differential cross section for the reaction 1.6 
Mev was obtained with targets and 45° the incident beam. 
The runs were interspaced with runs establish the tritium 
content the targets the same time and the same target position. The 
counter solid angle was 0.18 steradian measured from the target the 
the scintillation counter, 4.8 centimeters inside the front 
face the counter. This effective center distance was determined roughly 
calculating the distance inside the crystal which half the total absorption 
occurred and more accurately finding the intercept the distance axis 
for plot the reciprocal the square root the counting rate versus 
distance from the front end the crystal for number y-ray sources. 
The use center’’ permits one define efficiency factor 
independent the distance the source accuracy better than for 
distances the source about equal the length the crystal. The efficiency 
the counter was approximately for y-rays the energy range 
from 2.5 Mev, where this efficiency the ratio the number counts 
the spectrum above bias equal half the y-ray energy the number 
entering the counter solid angle. general, runs energies other than 
1.6 Mev were taken just prior just after runs 1.6 Mev. Runs taken 
with targets and agreed within 10%. Several runs taken target 
with sodium iodide counter in. long in. diameter agreed within 
15% with yield data from the larger counter. The results are collected 
Table III and Figs. and The first column the table lists the mean 
alpha-particle energy the target, the next three columns the 90° differential 
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Fic. The differential cross section measured 90° the incident beam. 
The solid vertical lines indicate the statistical uncertainties while the dotted lines indicate 
the estimated total standard errors. 
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Fic. center mass cross section factor for the reaction. Curve shows 
the present results and curve shows approximately the value obtained drawing 
smooth line through the points given Holmgren and Johnston (1959). 


TABLE III 
cross-section results 


total 

1.83 2.50 3.35+0.70 4.40+0.9 
1.32 1.89 2.74+0.50 3.58+0.6 
1.14 1.59+0.40 2.00+0.5 


cross sections, and the fourth column lists the total cross section integrated 
over all angles taking into account small correction for the nonisotropic 
yield higher energies. The final column lists the center mass value the 
cross-section factor used astrophysical work (Burbidge al. 1957) 


where the mean alpha-particle energy the target kev. The value 
Sasa function energy shown Fig. Within the accuracy constant 
with mean value 0.064 kev barns and mean square deviation 
from the mean. The estimated errors, however, are considerably greater than 
this and amount 25% shown Table III and Fig. These estimated 
errors are all cases greater than the statistical error about for points 
above 800 kev and about for points below 800 kev. The main contributions 
the estimated error are 15% for target tritium content, for statistics, 
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for scintillation counter efficiency, and for counter solid angle. The 
beam current measurements were made with current integrator 
accuracy 2%, and the targets were run +100 volts prevent secondary 
electron emission. 


DISCUSSION 


The total cross section obtained from the present results factor 
2.5 higher than that obtained Holmgren and Johnston (1959) over most 
the energy range measured, and there not much evidence for increase 
the factor the bombarding energy decreased. should pointed 
out that small errors estimating the mean energy the alpha particles 
the target cause increasingly large errors the value the energy 
lowered because the strong energy dependence the penetrability. Using 
our ratio and our total 0.064 kev barns, the factor for the transition 
the ground state 0.046 kev barns surprisingly good agree- 
ment with the value calculated Christy and Duck (1961), 0.05 kev 
barns. 

The observed ratio the number gamma rays the first excited state 
the number the ground state 0.45 1.6 Mev and 0.36 0.8 Mev 
good agreement with the value 0.42 obtained Tombrello and Phillips 
(1961) using direct capture theory. addition, the theory given these 
authors combined with the present cross section leads reduced width 
0.12 for the configuration alpha particle plus triton and, pointed 
out these authors, this much larger than the reduced width which 
supports the assumption two-particle cluster structure for Li’. detailed 
analysis the angular distribution data, particularly higher bombarding 
energies where there are pronounced departures from isotropy, have been 
undertaken yet. The main process consisting the capture S-wave alpha 
particles form relative p-state the alpha and triton with electric dipole 
emission leads isotropic distribution characteristic the results low 
bombarding energies. Tombrello and Phillips have calculated the contribution 
due magnetic dipole radiation following p-wave capture and find the 
the total cross section given Table III for 1.8 Mev. More detailed 
calculations are now being carried out check further into this matter. 
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THE EFFECTS ENERGETIC TRAPPED PARTICLES 
MAGNETOSPHERIC MOTIONS AND IONOSPHERIC CURRENTS! 


FEJER? 


ABSTRACT 


shown that the temporary distortion previously stable belt energetic 
charged particles, trapped the earth’s magnetic field, results system 
ionospheric currents and magnetospheric motions. The currents are caused 
the charge separation which occurs the course the subsequent adiabatic 
motion the energetic particles. The resulting current system closely resembles 
the well-known current system. The polarization fields, which drive these 
currents, give rise magnetospheric motions similar those Axford and 
Hines and the implications their work then also result. The distortion the 
belt could caused either directly, the introduction new group trapped 
particles mainly the day side the magnetosphere, indirectly, through 
the relatively sudden compression the earth’s field the day side, caused 
the arrival neutral stream charged particles from the sun. 


INTRODUCTION 


most recent theories magnetic storms the sequence observed events 
attributed the arrival neutral stream charged particles from the 
sun great intensification such stream. The changes the magnetic 
field, which take place during magnetic storm, are often interpreted terms 
uppér atmospheric current systems, although the location these external 
current systems cannot uniquely determined from measurements the 

Thus the decrease the field observed the ground during the main 
phase magnetic storm has been interpreted Schmidt (1924) and later 
Chapman and Ferraro (1931) ring current the equatorial plane 
distance few earth radii. Such ring current was postulated earlier 
Stormer (1911) explain the position the auroral zones with his theory. 

Likewise the diurnal changes the storm field, averaged over many 
magnetic storms, have been interpreted current system about 
height, the current system (Chapman 1935). Similar current systems were 
obtained for magnetic bays Silsbee and Vestine (1942). 

Singer (1957) suggested that the equatorial drift (westward for ions and 
eastward for electrons) energetic charged particles the cause the ring 
current and the main phase. The presence energetic charged particles has 
been subsequently confirmed observations with the aid satellites (Van 
Allen and Frank 1959). 

Dessler and Parker (1959) proposed that the trapping such particles 
made possible instabilities the boundary separating solar stream from 
the plasma confined the earth’s magnetic field. The instabilities make the 
penetration and subsequent diffusion neutral plasma blobs into the earth’s 
field possible relatively short time. 
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generally assumed that the plasma surrounding the earth confined 
the earth’s magnetic field distances the order earth radii. 
Measurements with the aid Explorer and Pioneer (Smith al. 1960) 
show that the earth’s field does extend distances this order but that there 
are systematic differences between the values the measured field and the 
field extrapolated from the ground. These differences were tentatively inter- 
preted Smith terms toroidal ring current about million 
amperes distance about earth radii from the center the earth 
the center the current. 

Gold (1959) has used the word ‘magnetosphere’ describe the plasma sur- 
rounding the earth and dominated the earth’s magnetic field. discussed 
the type motions possible the magnetosphere. 

Axford and Hines (1961) have examined the nature magnetospheric 
motions caused external neutral stream charged particles. They 
conclude that viscous-like transfer momentum from the external stream 
causes certain pattern motions inside the magnetosphere, and that this 
pattern significant many geophysical phenomena. They conclude 
particular that current system similar the system follows from these 
magnetospheric motions. 

the present paper mechanism described which leads ionospheric 
current systems similar the current system and pattern magneto- 
spheric motions similar that Axford and Hines (1961). The mechanism 
based the assumption that occasionally the axial symmetry the 
distribution energetic charged particles trapped the magnetosphere 
temporarily disturbed (for simplicity the earth’s unperturbed magnetic field 
regarded here axially symmetrical) and that symmetry re-established 
after some hours mainly the same equatorial drifts energetic charged 
particles (with spread velocities appropriate the energy spectrum) 
which according Singer (1957) causes the ring current and the main phase. 
Such periods temporary asymmetry may well occur several times suc- 
cession during magnetic storm. The ionospheric currents and magnetospheric 
motions occur only while the asymmetry persists, although ring current 
great distances continues flow after the re-establishment symmetry. 

Two possible causes temporary asymmetry are suggested the present 
paper. The first these the preferential trapping group energetic 
charged particles the day side. The second possible cause taken 
relatively sudden distortion the earth’s magnetic field solar stream. 
Although this latter process destroys the initial symmetry both the magnetic 
field and the particle distribution, its effects are shown similar those 
the first process which the symmetry the magnetic field least 
initially preserved. 

The nature possible motions the magnetosphere first discussed from 
general point view Section this paper. The consequences 
initial asymmetry the distribution trapped energetic charged particles 
with respect the magnetic field are examined Sections and with the 
aid the concepts developed Section 


FEJER: ENERGETIC TRAPPED PARTICLES 1411 


MOTIONS THE MAGNETOSPHERE 


There seem two different types ionized component gases present 
the magnetosphere addition certain amount neutral gas. The 
first type gas relatively low energy; its temperature believed 
about Its density known approximately distances 
about earth radii. The density below the peak the layer obtained 
from ionosonde data. heights between the layer peak and about 1000 
the density obtained from satellite observations different type and from 
observations incoherent scattering. still greater heights earth 
radii the observation nose whistlers (Pope 1961) yields electron densities. 

addition this low energy (thermal) component there are energetic 
trapped charged particles, some which have been observed with the aid 
satellite-borne counters (Van Allen and Frank 1959). Some these particles 
are believed solar origin; the consequences preferential trapping 
groups such particles the day side will discussed Section 

The motions the low energy component the absence energetic charged 
particles will considered first. follows from the frozen field concept that 
only such motions the low energy component the magnetosphere are 
permitted which ionized matter, which one time occupies tube force, 
continues all later times occupy tube force equal flux content. 
Gold (1959) has pointed out, the earth’s field must remain almost unaltered 
the motions. 

The motions the low energy component can have number possible 
causes. One these the tidal wind system the neutral air which believed 
cause the dynamo current system. Associated with the dynamo current 
system there electric field due polarization charges and derivable from 
potential. The global potential distribution over the thin spherical shell, 
which the dynamo currents flow, obtained from the dynamo theory (for 
example, Fejer 1953); may extended from this shell (the dynamo region) 
into the magnetosphere because the magnetic field lines are also lines 
equipotential.* knowledge the electric polarization field assumed 
then the drift velocity ionization across the lines force given 
where the magnetic induction the earth’s field. Motions 
ionization along the lines force are not considered this paper. 

addition the tidal wind system there probably less obvious wind 
system caused the rotation the neutral atmosphere with the earth 
dynamo region heights (about 90-150 km). The geometry this wind system 
would such that the electric field induced the moving air could balanced 
exactly polarization field and therefore ionospheric currents would 
flow. The polarization field would again cause drift velocity which 
would ensure the rigid rotation the magnetosphere with the earth. 


*It also follows from this condition that the potential distribution over the shell must 
symmetrical about the magnetic equator; the symmetry maintained with the aid magneto- 
spheric currents along the field lines between the northern and southern hemispheres. 
past treatments the dynamo theory these currents vanished account assumed 
symmetry about the magnetic equator. Similar assumptions symmetry are made the 
present paper. 
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follows from the type arguments used the dynamo theory that 
the neutral air the dynamo region were stationary then motions the 
low energy component across the lines force could occur. These arguments 
are only approximately correct since they neglect magnetospheric currents 
perpendicular the lines force. However, possible show that for the 
low pressures thermal ionization these perpendicular currents are very weak 
and the motions caused them are insignificant. 

The situation changes completely when energetic particles are present. 
Even the density these energetic particles were considerably lower than 
the density the low energy component, currents the order few million 
amperes could flow perpendicularly the field distance few earth 
radii. The presence such currents the basis Singer’s (1957) explanation 
the main phase. 

steady ring current caused energetic trapped particles, whose distri- 
bution symmetrical about the earth’s magnetic axis, free divergence and 
therefore does not lead ionospheric currents. The next section will show 
how lack such symmetry leads ionospheric currents and magnetospheric 
motions. 


THE EFFECTS ASYMMETRIC TRAPPING PARTICLES 


convenient start with the consideration relatively simple type 
asymmetry which new group trapped particles introduced pre- 
dominantly the day side the magnetosphere. For simplicity the resulting 
ring current assumed line current the equatorial plane. The value 
the current taken greater the day side than the night side 
with abrupt changes current the plane. 

the subsequent discussions the ring current taken caused entirely 
particle drift. There are additional currents due the gyration particles 
but they will ignored because they are free divergence. 

The use line current instead current with finite and large cross 
section, and the assumption abrupt rather than gradual changes, are 
idealizations which simplify the illustration the mechanism without changing 
its essential nature. The assumption initial trapping the day side not 
justified here theoretically; intuitively appears reasonable. The same 
assumption has been made Chamberlain (1960), who, however, regard 
subsequent departures from neutrality unimportant. Such departures from 
neutrality are essential the present mechanism. 

Since axially symmetrical ring current does not give rise ionospheric 
currents, sufficient consider only the excess ring current for the 
calculation the initial ionospheric current system. This excess current 
flows along half circle the equatorial plane the day side, from sunset 
sunrise, illustrated Fig. the meridian plane the 
low energy component forms highly conducting paths along the lines force 
which connect the equatorial current with the polar regions. The circuit 
then completed two ionospheric currents the dynamo region, 
which flow from the northern and from the southern 


e 
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TO SUN 


Fic. view magnetospheric and ionospheric currents associated with 
excess ring current the day side. 


hemisphere. These ionospheric currents are schematically indicated single 
flow lines which pass through the north pole and the south pole re- 
spectively. 

The divergence the current system flowing the ionospheric shell 
thus known. the idealized current system Fig. the divergences are 
points and the northern and and the southern hemisphere. The 
ionospheric currents are driven quasi-static electric field which has polari- 
zation charges its sources and which derivable from potential. The 
determination such current system known divergence, driven 
irrotational electric field, problem well known from the dynamo theory 
where has been attacked numerical methods (for example, Fejer 1953, 
Maeda 1955). 

With the simple geometry Fig. and assumed that the conductivity 
the ionospheric shell does not vary much with geographical position over 
the polar regions, the general nature the resulting ionospheric current 
system may obtained from relatively simple considerations. 

There are two types current the ionospheric shell. The first type 
parallel the horizontal component the electric field while the second type, 
the Hall current, perpendicular it. the conductivity the ionospheric 
shell approximately independent position then the Hall current almost 
free divergence, and the divergence the parallel current proportional 
the divergence the horizontal component the electric field. Thus the 
horizontal component the electric field will have positive divergence, 
source and sink the northern hemisphere; the southern 
hemisphere the source and the sink. The Hall current caused 
such electric field flows loops around the points and 
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indicated Fig. account the nature ionospheric conductivity (for 
example, Baker and Martyn 1953) the Hall current considerably greater 
than the current which flows parallel the electric field and completes 
the circuit the magnetospheric currents. (The ionospheric current flows 
directions perpendicular the Hall current loops. For simplicity only 
single flow line this current indicated for each hemisphere; likewise only 
one the many Hall current loops around each the points and 
shown.) 

The Hall current loops Fig. show some resemblance the current 
system Chapman (1935) and the current systems derived from magnetic 
bays Silsbee and Vestine (1942). 

The ionospheric currents described above are driven polarization field 
The same polarization field also causes motions the magnetosphere 
velocity EXB/B?. The Hall current results from the motion the electrons 
alone this velocity the dynamo region, where the motion ions stopped 
collisions while the electrons are still free drift. picture the motions 
ionization above the dynamo region thus obtained simply reversing 
the arrows the current loops Fig. The magnetic field lines extend the 
motions into the magnetosphere the previously described manner. 

The system currents and motions Fig. only represents the initial 
effect the excess energetic particles the day side. The energetic particles, 
both newly introduced and previously present, are themselves affected the 
electric fields which cause the magnetospheric motions. The effect the 
previously present belt energetic trapped particles particularly important 
since the ring current carried likely very much greater than the 
the center the belt towards the night side the magnetic north pole 
shown Fig. Figure view the earth from above the north 
magnetic pole. The position the belt trapped particles indicated pro- 
jection onto the northern hemisphere along the lines force. The shaded area 
between the two concentric circles represents the displaced belt. assumed 
for simplicity that such displacement occurs first and that then followed 
the usual equatorial drift the particles the belt; reality the two 
processes occur simultaneously. 

The arguments leading the currents and motions Fig. can now 
used again for the rather different type asymmetry represented Fig. 
the magnetosphere drifts the ring current type, indicated the lines 
with arrows (which are arcs circles with their center marked for 
protons and for electrons, are not free divergence. Currents flow down 
the lines force the dynamo region along the edges the belt marked 
with the signs. The currents then return along other lines force along the 
edges marked with the signs. previous arguments the same and 
signs represent the divergence the horizontal component the electric 
polarization field the ionospheric shell. The currents parallel the electric 
field flow from the the signs while the Hall currents flow along the 
lines equipotential the shell. Only the Hall currents, which are more 
intense than the parallel currents, are indicated schematically Fig. 
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SUN 


Fic. The ionospheric Hall currents (interrupted lines) caused temporary displace- 
ment the trapped particle belt, seen from above the magnetic north pole The pro- 
jection the belt along the magnetic lines force onto the northern hemisphere indicated 
the shaded area between the two concentric circles. The center these circles dis- 
placed from towards the night side. 


noted that the electric field, determined the divergences shown 
Fig. the and signs, must particularly intense narrow 
ionospheric zone which obtained projection the particle belt along the 
field lines. Intense Hall currents must flow the same zone. 

There will again motions ionization the magnetosphere the 
velocity picture the motions immediately above the dynamo 
region may again obtained reversing the arrows the Hall current 
loops. The velocity the motions will greatest above the same zones where 
the Hall currents are most intense. 

The current system Fig. resembles the current system even more 
closely than the current system Fig. because predicts intense currents 
relatively narrow zones (the auroral electrojets) and because predicts 
current loops different sense for high and low latitudes. 

consideration Fig. should make clear that the previous presence 
several narrow belts trapped particles would lead considerably more 
complex current system, particularly near the projection the belts where 
complete reversals the direction the currents could occur within relatively 
short distances. 


The interpretation given the preceding discussion the current system 
Fig. implies that the resultant current system obtained some sort 
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superposition Fig. and Fig. The arguments used the derivation 
these current systems were semiquantitative nature. more quantitative 
theory would have make more precise assumptions about the energy 
spectrum and spatial distribution the previously present and the newly 
introduced energetic trapped particles well about ionospheric electron 
densities and collision frequencies. The theory would have describe the 
subsequent motion both the thermal and the energetic components the 
magnetosphere until equilibrium again established largely the ring 
current type drift energetic particles with spread drift velocities 
appropriate their energy spectrum. 

the energy the electrons were predominant then eastward rotation 
the initial current system could expected with the subsequent weakening 
the currents caused the spread drift velocities. 

the energy the positively charged particles were predominant then the 
rotation would westward observed from co-ordinate system which 
rotates with the earth. If, however, co-ordinate system, whose orientation 
remains fixed with respect the sun, were used then the rotation the 
current system could either eastward westward depending whether 
the eastward velocity produced the earth’s rotation (by means the 
electric polarization field discussed Section greater smaller than the 
ring current type westward drift velocity. the ring current type westward 
drift velocities the positive energetic particles were about equal magnitude 
their eastward velocity caused the earth’s rotation then the current 
pattern would hardly rotate all and would merely decay account the 
spread drift velocities. 


THE EFFECTS ASYMMETRIC CHANGE THE 
MAGNETIC FIELD 

asymmetry similar that Fig. could produced entirely 
different manner through relatively sudden compression the earth’s 
magnetic field the day side solar stream originally suggested 
Chapman and Ferraro (1931) their explanation the initial phase 
magnetic storm. The compression assumed occur time much longer 
than the gyration periods the particles but much shorter than the time 
required the particles once around the earth their ring current type 
drift motion. 

The result such compression the field has been derived Parker 
(1960) whose Fig. resembles the present Fig. the asymmetry Fig. 
were produced the manner suggested Parker, then the mechanism the 
present paper would not require the introduction any additional trapped 
particles. The resulting initial current system would then that Fig. 
alone. 


GEOPHYSICAL CONSEQUENCES MAGNETOSPHERIC MOTIONS 


The motions following from the mechanism the present paper are similar 
those Axford and Hines (1961). Their discussion some the conse- 
quences these magnetospheric motions will not repeated here detail. 
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The motions observed visual and radio aurorae are perhaps the most direct 
these consequences. 

more indirect effect the motions the acceleration process ‘adiabatic 
compression’, which occurs the night side and was put forward Axford 
and Hines possible basis for auroral theory. equivalent explanation 
this acceleration process has been given their paper terms the 
motion individual energetic charged particles. 

This acceleration process the night side with the resulting lowering 
mirror points and eventual loss trapped energetic protons and electrons 
(at different latitudes) could lead additional current system. The forma- 
tion current system similar manner, but different acceleration 
process, has been suggested Chamberlain al. (1960) and Kern (1961). 


CONCLUSIONS 


was attempted show the present paper that the assumption 
certain types temporary asymmetry the distribution trapped particles 
relative the magnetic field leads system currents, motions, and 
acceleration processes which appear explain, least qualitatively, many 
the observed phenomena associated with magnetic disturbances and 
aurorae. 
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TRANSITION RANGES FOR THREE-DIMENSIONAL 


MAGARVEY AND Roy BISHOP 


ABSTRACT 


The motions set liquid due the passage drop immiscible 
liquid have been examined for number liquid-liquid systems. Wake con- 
figurations corresponding Reynolds numbers between and 2500 have been 
photographed, and the transition ranges between the various wake configurations 
have been determined. Several systems were used obtaining data 
wake patterns, and was found that identical patterns correspond approxi- 
mately the same Reynolds numbers regardless the system. However, the 
exact transition number cannot determined for the general case drop 
deformations depend the physical characteristics particular system. 
The range Reynolds numbers corresponding transition from one wake 
pattern another somewhat less than 20. 


INTRODUCTION 


Apart from practical considerations, wakes systems are 
interesting physical phenomena their own right. The periodicity and 
symmetry evident some the three-dimensional wake forms are com- 
parable with any associated with two-dimensional patterns. Three-dimensional 
wakes corresponding drop motions characterized Reynolds numbers less 
than about 210 are single threads, and those corresponding numbers 
between 290 and 700 consist double row vortex rings. This double row 
configuration has been photographed Magarvey and Bishop (1960, 1961). 

The Reynolds number the most convenient flow index which wake 
configurations associated with systems can compared, but 
due the deformable nature the drop the precision associated with solid— 
liquid systems lacking. The Reynolds number given 


(1) 


where some characteristic length, the relative velocity between the 
drop and the quiescent phase, and the kinematic viscosity the disperse 
liquid. these experiments, the major axis the deformed drop was found 
more pertinent dimension than the equivalent diameter. 

two-dimensional wakes have been investigated for wide range 
Reynolds numbers, brief examination the phenomena indigenous 
consideration the three-dimensional counterparts. Essentially, with increas- 
ing Reynolds numbers, increasing amount vorticity originating the 
boundary layer concentrated into two eddies the immediate downstream 
side the obstacle. The vortices thus formed gain strength, increase size, 
and change shape the Reynolds number increased. Eventually the 
standing eddies become elongated that they break down, and, the case 
cylinder, characteristic state flow developed which vortices are 
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shed alternately from the two sides the obstacle. The eddies are shed 
regular intervals, forming wake configuration called vortex street. 
general, the Reynolds number which the shedding initiated depends 
the shape the obstacle. The frequency shedding given terms 
dimensionless parameter called the Strouhal number. This number given 


(2) 


where the frequency, the velocity, and the same dimension 
equation (1). Rayleigh (1915) obtained empirical relationship between the 
Strouhal and Reynolds numbers for flow past circular cylinders. The equation 
relating these two dimensionless groups and the range for which they are 
valid are indicated 


The theoretical three-dimensional generalization the standing eddy 
phenomenon ring with the axis coinciding with the direction flow. This 
ring has been observed the downstream side flat plates and bodies 
revolution for limited range Reynolds numbers. Stanton and Marshall 
(1930) observed vortex rings for very small Reynolds numbers behind circular 
plates. The nature the wake outside the vortex sheet enclosing the two- 
dimensional eddy pair the three-dimensional ring has received little attention 
other than for cases periodic shedding. This report concerns three-dimen- 
sional wakes corresponding range Reynolds numbers extending from 
about 2500. The range numbers over which the transition from one 
wake pattern another occurs has been found experimentally less 
than 20. 


WAKE CLASSIFICATION 

convenient classify wakes according the nature the motions 
set the continuous phase far behind the moving body. Although the 
motions within the vortex sheet the downstream side the body determine 
the nature the trail far behind, these motions are much more difficult 
observe. systems, six distinct and reproducible trail configura- 
tions are recognizable for the entire range stable drops. The range 
numbers corresponding each classification indicated Table The line 


TABLE 
Range 
Class Reynolds numbers Nature trail 
0-210 Single thread 

210-270 Double thread 

270-290 Double thread with waves 
290-410 Procession vortex loops 

290-700 Double row vortex rings 


700-2500 Asymmetrical wake 
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demarcation not sharp between any two classes the values are result 
observations different systems. The identical lower limit for classes 
and cannot explained from examination the photographic evidence. 


STEADY WAKES 

Two flow types are associated with the single thread wake. For very small 
Reynolds numbers the streamlines follow the interface and exhibit symmetry 
for spheres and other symmetrical forms. For Reynolds numbers above 
critical value, the streamlines separate from the surface and form closed 
region immediately behind the moving body. The exact transition Reynolds 
number has not been determined, but experimental evidence indicates value 
less than 10. The trail left the passing drop single thread regardless 
separation, provided the Reynolds number less than 210. The vortex layer 
which leaves the interface forms vortex sheet separating the wake liquid 
from the main flow. Inside the sheet inflow along the axis the wake gives 
rise vortex ring immediately behind the drop. The single thread sufficient 
maintain equilibrium between the rate which vorticity generated 
and the rate which dispensed into the main stream. Figures and 
show the single thread some distance behind the advancing drop and section 
the ring inside the sheet. 

One the most interesting wake configurations the double thread observed 
behind drops characterized Reynolds numbers between 210 and 270. 
asymmetry the circulatory motion within the vortex sheet produces 
corresponding asymmetry the sheet itself and resultant deviation from 
vertical line fall. rotation observed, assumed that the lack 
symmetry responsible for the sidewise force component. The exact nature 
and cause this lift cannot inferred from the photographic data. The 
double thread maintains equilibrium between rates generation and 
diffusion vorticity. 

Gunn (1949) observed deviations from the vertical direction fall for 
water drops the same Reynolds numbers characterize these class 
drops systems. The water drops were falling through air, and 
Gunn attributed the behavior mechanical resonance phenomenon 
associated with eddy detachment and natural frequency drop oscillation. 
Since comparable motion can observed directly systems, 
entirely possible that the two phenomena are identical. However, there 
are definitely aeolian sheddings and oscillations associated with liquid— 
liquid systems this range Reynolds numbers. Figures and indicate 
the wakes far behind and close the moving drop. The circulations within 
the vortex sheet are evidence Fig. The transition Reynolds number 
from single double thread configuration fairly sharp for specific 
system, but for the general case, 210 the approximate value. 


PERIODIC WAKES 
The state equilibrium set under conditions corresponding the first 
two classes disturbed increasing further the Reynolds number. The rate 
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behind carbon tetrachloride drop moving through water 


terminal velocity. Reynolds number 170. 
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(a) and (6) Double thread wakes for two drops. The dye concentration 
the pattern circulation within the vortex sheet. Reynolds numbers 225 and 220 


respec live 


. 
& 
f 
j 


(a) 


and (b) Wavy double thread wake behind drop monochlorobenzene moving 
through water. 
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(a) and (b) Two views class wake. Interval seconds between 
photographs. Reynolds number 370. 
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(a) (b) 


(a) and (b) Two views single cycle loop wake. Circulation inside vortex 
sheet shown from two angles. 
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Class wake consisting double row vortex Reynolds number 
Kings associated with class wake. Reynolds number 
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which vorticity diffused from the sheet into the main body the fluid 
remains practically constant, but the increased rate which transferred 
the vortex ring creates unstable conditions within the vortex sheet. The 
release vorticity regular intervals the first suggestion periodic wakes. 
Basically, the process one build-up and release, but sizeable portion 
the ring escapes through opening the end the vortex sheet during the 
cycle. When the vortex strength the ring reaches critical value, sudden 
motion the ring disturbs the sheet, which, turn, responsible for release 
vorticity and consequent return the ring its original position and 
shape. The class III wake corresponds narrower range Reynolds 
numbers than the first two, and could considered transition range. 
However, according the system classification based the nature the 
trail far behind the drop, this type wake must considered separate 
class. The wake appears double thread characterized equally spaced, 
wavy disturbances. Figures and indicate the wavy trail and the 
circulation pattern immediately behind the drop. 

For Reynolds numbers greater than 290, the trail assumes repetitive 
pattern and the circulation within the sheet ceases resemble symmetrical 
vortex ring. Two distinct wake configurations correspond the range 
Reynolds numbers extending from 290 410: the first characterized 
series vortex loops connected intricate system vortex filaments, 
and the second the double row vortex rings described Magarvey and 
Bishop. The two wake types have common lower limit, but the symmetrical 
vortex ring pattern extends upper limit 700. The occurrence either 
the loop the ring wake the overlapping region Reynolds numbers 
appears matter chance. Initial conditions drop formation probably 
are responsible for the flow which determines the nature the wake. 

the cycle build-up and release, the vorticity generated the boundary 
layer concentrated diametrically opposite sides the flow axis within 
the vortex sheet. The sections which the vortex strength greatest are 
alternately discharged into the main body the fluid. With each ejection 
portion the sheet carried away. The vortex element discharged into the 
stream interacts with the disperse liquid form regular wake pattern. The 
vortex loop configuration termed class wake and the vortex ring street 
class wake. Figures and show two separate loop wakes photo- 
graphed different angles. Figures and (6) show the circulation patterns 
within the vortex sheet for this type wake. 

examination many loop wakes indicates that the frequency discharge 
function velocity and the physical characteristics the disperse phase. 
The physical properties the drop phase are less importance, but they 
appear factors determining the velocity the drop. Figure shows the 
Strouhal number plotted against the Reynolds number for wakes over the 
entire range Reynolds numbers for which class wakes are observed. 
the best fitting curve straight line, the relationship between the 
and velocity given terms the slope. 


t 
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where the kinematic viscosity the field phase. Magarvey and Bishop 
(1961) found linear relationship between the Strouhal and Reynolds numbers 
for the range corresponding class wakes. (1938) made comparable 
plot for rigid spheres covering range Reynolds numbers for which class 
wakes are observed systems. However, the wakes behind the 
spheres investigated not exhibit the periodicity the class 
wakes behind drops. 

Figure shows class wake photographed carbon 
water system. The frequency discharge greater than those class IV. 

The symmetry and periodicity characterizing the wakes classes and 
are not evidence for Reynolds numbers greater than 700. Vortex rings are 
still the main features these wakes, but there appears preferred 
plane which the individual elements are discharged into the quiescent phase. 
The vorticity generated the boundary layer concentrated into the liquid 
within the sheet such rate that definite pattern circulation can 
established before element discharged. The discharge one element does 
not initiate conditions which lead discharge from the diametrically 
opposite side the wake axis. Figure shows typical class wake. 


0.09 


0.07 


0.05 


STROUHAL NUMBER 


280 300 320 340 360 380 400 


Fic. showing relationship Strouhal and Reynolds numbers for loop wakes 
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THE STATISTICAL PROPERTIES 
OPTICAL AND RADIO AURORA! 


ABSTRACT 


The latitudinal variation certain auroral forms has been investigated 
means all-sky camera records made central Canada during the The 
forms chosen were those that may have significance theories auroral morpho- 
logy. The occurrence seems concentrated near the auroral 
zone. One type radio reflection appears associated with the curl 
forms and another type with patchy, early morning aurora, which occurs 
most frequently well south the auroral zone. 


INTRODUCTION 

The remarkable similarity between certain auroral configurations and the 
instability patterns that develop sheet beam charged particles moving 
along magnetic field, seems first have been brought out discussions 
between Webster and Kavadas (see Webster 1957). The auroral configuration 
develops kink which rapidly progresses into curve (as seen plan 
view) and thence into nearly isolated spiral configuration. The curl form 
appears represent the transition between the quiet arc form and the broken 
ray forms. the obvious similarity between this auroral sequence and the 
break-up sheet beams more than superficial, implies that the auroral 
configurations are actually determined beams incoming particles that 
travel considerable distances. The alternative view that the shapes auroral 
forms are determined local conditions the ionosphere. 

Recently, Kern and Vestine (1961) have examined several features auroral 
morphology, including motions auroral forms, the assumption that 
sheet beams are involved and that the sheet beams are made particles 
trapped the geomagnetic field. Earlier, Gartlein (1959) pointed out that 
the auroral configurations are governed the precipitation geomagnetically 
trapped particles, the occurrence curl forms should vary with geomagnetic 
latitude, because the varying length the geomagnetic field lines. 
might expected that the relative occurrence curl forms should increase 
with increasing geomagnetic latitude, least some latitude where the 
auroral arcs cease occur. 

the course radar investigation aurora, Lyon (1960) found notable 
similarity between the diurnal variation curl forms and the occurrence 
radar echoes for particular location (near Flin Flon). discussing this 
result, Forsyth (1960) drew attention the particular geometry involved 
which made uncertain whether the association curls with radar echoes 
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was due increase ionization when the curls are formed kind 
observational selection that was peculiar the radar situation. was sug- 
gested that this point could resolved means continuous-wave, bistatic, 
radio observations the same general region. 

This paper attempts shed some light these two characteristics curl 
forms: the variation occurrence with latitude and the association with 
radio reflections. the course the study interesting characteristics were 
found also for the patchy, early-morning aurora. The results were obtained 
analysis all-sky camera films and bistatic radio records obtained during 
the International Geophysical Year. 


RESULTS AND ANALYSIS 

During the all-sky auroral cameras were operated Uranium City 
(59.5° N., 108.6° W.), Flin Flon (54.8° N., 101.8° W.), Saskatoon 
106.6° W.), and Regina (50.4° N., 104.7° W.). The choice suitable period 
for analysis required compromise between the need for statistically 
acceptable body data and the requirement that should homogeneous 
possible. The vagaries weather and equipment performance precluded 
any completely satisfactory solution. The period chosen was November, 
1957, February, 1958. The number useful records obtained this 
period for each station shown Table The number hours clear- 


TABLE 
All-sky film records 


No. hours Percentage 


useful film dark hours 
Uranium City 183 
Flin Flon 437 
Saskatoon 286 
Regina 368 


weather film obtained, particularly Uranium City, was small fraction 
the total number dark hours. However, this situation prevailed throughout 
the I.G.Y. and there was, fact, strong concentration useful records 
the period chosen. felt that, while the number data less than 
desirable, the validity the analysis not likely improved the 
inclusion additional, much less homogeneous, data from the remainder 
the 

The films were analyzed assigning each 10-minute period classification 
representing the most prominent type aurora present. The classifications 
used were: 

(i) arcs, including only well formed arcs, each with clearly defined lower 
edge without scallops folds; 

(ii) curls, the curl forms described earlier including arcs with wavy 
structure along the lower edges seen the photographs; 

(iii) patchy, including all the forms that appear luminous patches the 
photographs but probably referring mostly pulsating patches; 
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(iv) other, all other forms including diffuse surfaces, active rayed bands, 


and rays. 


The results the analysis are presented Fig. where the percentage 
occurrence (percentage the total number clear-sky, 10-minute periods) 


(a) ARCS 
URANIUM CITY 
FLIN FLON 
SASKATOON 
REGINA 


N 
o 


(c) PATCHY AURORA 


OCCURRENCE (% TIME) 
@ 


(b) CURLS 


UNIVERSAL TIME 


UNIVERSAL TIME 


Fic. The variation occurrence with time various auroral forms, for the dark hours. 
For each station, local midnight occurs near hours U.T. (see Table 


each classification each hour the night plotted. Except for the first 
and last hours U.T. and U.T.) there was fairly uniform distri- 
bution useful film records over the night hours for each station. 


Bistatic Radio Observations 


The radio measurements made use two the bistatic paths, 
Saskatoon and Yellowknife-The Pas, that are described elsewhere (Forsyth, 
Green, and Mah 1960). may seen from Table II, the geomagnetic latitude 


TABLE 


Geomagnetic latitude and local time 


Geomagnetic 


latitude 
Uranium City 67.4° 
Flin Flon 
Saskatoon 60.5° 
Regina 
Yellowknife-The Pas (mid-point) 


(mid-point) 64.9° 


Time local 
midnight (U.T.) 


min 
min 
6min 
min 
min 


the mid-point the path slightly greater than 
that Flin Flon. The mid-point the Yellowknife-The Pas path near, 
but south of, Uranium City. The records were the same those used earlier, 
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but this time they were subjected more detailed analysis. before, the 
record interval used for analysis corresponded minutes time. For each 
interval the presence enhanced signal, due scattering the very high 
frequency radio waves auroral ionization was listed event. The earlier 
study made use only and events. The present study included also the 
events classified and The classifications used here and A3) 
were described and illustrated earlier paper (Collins and Forsyth 1959). 
Although the events were low amplitude the records was found 
that with care they could distinguished with reasonable certainty. Since 
the analysis was confined the night hours, events were not encountered. 
indicated both the earlier studies, the occurrence events these 
latitudes negligibly small and will not considered further. 

The radio records were independent weather conditions and, except for 
occasional periods equipment failure, were continuous throughout the 
3-month period. all, the number hours for which useful records were 
obtained amounted 1200 for the path and 1233 for 
the Yellowknife-The Pas path. The results the analysis are presented 
the next section. 

DISCUSSION 

unfortunate that the all-sky camera results from Uranium City are 
based few data. The very high percentage occurrence curl forms would 
seem offer strong support Gartlein’s contention that the occurrence 
these forms should increase with increasing geomagnetic latitude. While the 
details the hour-to-hour variations for Uranium City, shown Fig. 
should treated with some caution seems clear that curl forms did occur 
City substantially greater numbers than Flin Flon. There 
was also tendency for these curl forms concentrated near the middle 
the night. 

contrast those from Uranium City, the films from Flin Flon and 
Regina were excellent and those from Saskatoon were acceptable. The graphs 
Fig. (6), and (d) show, quite clearly, marked decrease occurrence 
aurora with decreasing latitude. The percentage occurrence all aurora 
each station was: Uranium City, 88%; Flin Flon, 65%; Saskatoon, 52%; 
and Regina, 30%. Even more interesting the sharp relative decrease 
occurrence the curl forms. would appear that the curl forms only occur 
significant numbers near the auroral zone, conclusion that accord 
with the implications the sheet-beam mechanism curl formation discussed 
Gartlein. 

Another interesting feature Fig. the extremely frequent occurrence 
the aurora the morning hours. While the occurrence this type 
aurora has often been noticed visual observers, seems have received 
little attention. Heppner (1954) showed that was primarily early-morning 
phenomenon and that seemed occur most frequently near geomagnetic 
latitude 64° (well south the maximum the auroral zone). Both these 
conclusions are supported the present results. 
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unexpected feature Fig. the relative unimportance the arcs 
(Fig. This result not agreement with that Lyon, who found 
high peak the occurrence quiet arcs during the early part the night. 
The difference between the two results probably due the more stringent 
criteria used the present study. examination Fig. does show 
tendency for the arcs occur later times with decreasing latitude. This 
characteristic also supports the sheet-beam hypothesis. 

The results the analysis radio records are shown Fig. Also shown 
the same figure are the diurnal variations curls and patchy aurora 
Flin Flon and Uranium City (from Fig. and (c)). spite the separations 


CURLS RADIO-AURORA 
URANIUM CITY o--o Y'KNIFE-THE 


FLIN FLON 


PATCHY AURORA 

URANIUM CITY 
RADIO-AURORA 


Fic. comparison the variation occurrence with time curl forms and one 
type radio scattering and similar comparison for patchy aurora and another type 
radio scattering 


between the all-sky cameras and the path mid-points and the uncertainties 
the results from Uranium City, seems clear that there close similarity 
between the diurnal variation the curl forms and that the events. 
This confirms the finding Lyon and indicates that his result was not due 
the particular geometry his experiment. must concluded that the 
observed association between the radio reflections and the curl forms due 
the presence intense ionization these structures. Fig. the radio 
events observed the path appear lag slightly 
time behind the curl forms Flin Flon. This lag probably due the fact 
that local time the radio path mid-point about half hour behind local 
time Flin Flon (see Table 
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Fig. there also some similarity between the variations the 
events and the patchy aurora. The association has often been noticed but this 
seems the first time that statistically significant evidence has been 
presented. The radio equipment used during the I.G.Y. was considerably less 
sensitive than that used for other investigations. When more sensitive equip- 
ment used, the events usually dominate the results. this had been the 
case for the I.G.Y. experiment the similarity between the diurnal variation 
the patchy aurora and that the events would have been even closer. 
Since the faint, patchy aurora that occurs the end optical display 
often disregarded visual studies, not surprising that many observers 
have concluded that there little relation between optical and radio aurora. 

the association between the curl forms and the radio reflections can 
accepted, then the radio evidence can used expand our knowledge 
auroral morphology. There ample evidence (Forsyth, Green, and Mah 
1960; Green 1961) show that the peak occurrence events occurs near 
midnight subauroral latitudes and progressively earlier with increasing 
latitude. From the discussion Kern and Vestine would appear that this 
behavior might expected the curl forms, they arise result 
instabilities set sheet-like configurations particles trapped the 
geomagnetic field. 
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THE ELECTRICAL RESISTIVITY Mn;ZnC 
BETWEEN 4.2° AND 630° 


ABSTRACT 


The electrical resistivity alloy having approximately the composition 
Mn;ZnC has been measured between 4.2° and 630° From high residual 
value the resistivity increases with temperature complicated way 
550° becoming temperature independent above that point. Two anomalous 
regions are evident corresponding the magnetic transitions occurring this 
substance ~240° and ~450° The general features the results are 
interpreted the basis spin-disorder scattering picture. 


INTRODUCTION 


Alloys with compositions near Mn;ZnC exhibit unusual magnetic properties 
(Butters and Myers 1955). Below Curie temperature (~450° K), which 
varies with exact composition, these substances become ferromagnetic. 
lower temperature, ~240° neutron diffraction (Brockhouse and Myers 
1957) and heat capacity (Swanson) data indicate another co-operative 
transition, the saturation moment this point passing continuously through 
maximum value. the basis their neutron diffraction studies, Brockhouse 
and Myers have concluded that between and all atoms these 
alloys are characterized the same moment and that these are ferro- 
magnetically aligned. Below the observations are described assuming 
that the atoms have different moments, Bohr magnetons, and 
that ferrimagnetic arrangement spins exists. 

view the redistribution spins among atoms presumed accom- 
pany the ordering transition have felt some interest 
examine the temperature variation electrical resistance over 
extended range including both and Such observations might 
expected suggest what role, any, the electrons responsible for the magne- 
tization play the conduction process. 


SPECIMEN MATERIAL 


alloy having composition close was prepared sintering 
and for days 600° and then for days 900° evacuated 
quartz capsule. The high temperature sinter yielded brittle alloy metallic 
appearance. diffraction photographs showed the alloy con- 
tain approximately 95% the desired perovskite structure with lattice 
parameter 3.921; The density the alloy was compared 
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with the density 6.63 calculated from the lattice parameter. According 
Butters and Myers (1955) such lattice parameter should characterize 
material the composition with Curie temperature 
416° Chemical analysis indicated, however, more nearly 62.8 at. man- 
ganese, suggesting that the average lattice parameter the specimen was 
lower than that quoted above and that the observed Curie temperature should 
higher than 416° crude determination the Curie temperature 
indicated, fact, the value ~453° 


EXPERIMENTAL METHOD 


The electrical resistivity was measured with accuracy 
0.5% the potentiometric method from 630° Mechanical and 
soldered contacts gave consistent results. Above the liquid nitrogen range, 
temperatures were measured +1° with calibrated 
thermocouple. Lower temperatures were determined means the vapor 
pressure temperature relation for the particular refrigerant employed, for 
example, liquid He. 


RESULTS AND DISCUSSION 
temperature Fig. The magnitude the resistivity 298° 370 
micro-ohm cm, roughly one-half that reported Butters and Myers for 
loosely sintered alloy similar composition. When made for the 
fact that the present material, while closer density the crystal value, 
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still not fully compact, its room temperature resistivity comparable with 
those other ferromagnetic metals and alloys. The general downward trend 
further indication the metallic character 

Several features the vs. curve are noteworthy. First all, becomes 
essentially temperature independent above about 550° Secondly, distinct 
anomalies appear the curve corresponding the two magnetic transitions 
The higher these anomalies temperature includes small 
but reproducible region negative slope. plot the temperature coefficient 
vs. shown Fig. brings out more clearly the features the 
anomalies. The maximum ~241°K coincides quite closely with the 
transition temperature determined calorimetrically Swanson similar 


160 240 320 400 480 560 640 
TEMPERATURE (°K) 


for Mn;ZnC. 


alloy, 238° The negative minimum near 450° occurs approximately 
the Curie point, estimated for this alloy from rough magnetic observa- 
tions. may remarked that the range which this material exhibits 
negative temperature coefficient resistivity much more restricted than 
that found Butters and Myers for their loosely sintered alloy. 

The general character the vs. curve for Mn3ZnC just described 
suggests that important resistive mechanism this substance probably 
spin-disorder scattering.* This scattering arises from the interaction the 
conduction electrons with unpaired electrons assumed localized the 
magnetic ions and largely responsible for the magnetization the 
solid. the ionic moments are presumed perfectly ordered 
that scattering occurs. Disorder among the ionic moments increases with 
rising temperature, producing corresponding increase the spin-disorder 
contribution the resistance until the transition temperature reached. 
Above the transition temperature long-range order completely destroyed 


review this subject with full references has been given Coles (1958). 
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and the absence short-range order the spin-disorder resistance would 
expected temperature independent. DeGennes and Friedel (1958) have 
shown that short-range effects can cause the resistivity not only tem- 
perature dependent immediately above the transition point but actually 
exhibit negative temperature coefficient and some cases pass through 
shallow minimum before achieving constant value. 

Two ordering transitions occur The knee the vs. curve 
near 200° (see Fig. suggests approach constant resistivity (~320 
micro-ohm cm) accompanying the lower transition. Complete spin disorder 
not achieved (~240° K), however, and the increase resistivity 
resumed the temperature raised through this point. Well above the 
upper transition temperature, the resistivity finally becomes constant about 
410 micro-ohm cm, but not before has passed through shallow minimum 
near tempting attribute this minimum the short-range order 
effects referred This possibility given added weight the fact 
that the specific heat anomaly observed Swanson actually exhibits 
very pronounced for usually clear indication the persistence 
short-range order. Interestingly enough, the specific heat anomaly 
has only small tail. This anomaly, however, larger overall than the upper 
one. While has not been possible separate completely these anomalies 
from one another from the lattice heat capacity, evident qualitatively 
that the gain both magnetic enthalpy and entropy, i.e. spin disorder, 
the lower transition point are larger than the corresponding quantities 
between and This consistent with the observation that between 
and 240° the resistivity increases larger amount than does between 
240° and 450° and further supports the assumption that spin disorder 
provides the dominant scattering mechanism this substance. 

Measurements the electrical resistivity between 4.2° and 
630° are thus apparently interpretable the basis the spin-disorder 
scattering picture. This suggests that the electrons whose spins are primarily 
responsible for the magnetization are localized manganese ions and that 
another group electrons whose number independent temperature carries 
the electric current. This appears true both above and below the spin 
rearrangement transition temperature There indication the present 
results that delocalization magnetic electrons occurs perhaps 
happens magnetite (Verwey, Haayman, and Romeijn 1947). 


*Dysprosium also exhibits two magnetic transitions with corresponding resistivity anomalies 


(Coles 1958). 
example such minimum that found for alpha-manganese its Néel point 
(Coles 1958; White and Woods 1957). 
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UNIFYING THEORY HIGH-LATITUDE GEOPHYSICAL 
PHENOMENA AND GEOMAGNETIC STORMS! 


ABSTRACT 


This paper concerned with the occurrence high latitudes large number 
geophysical phenomena, including geomagnetic agitation and bay disturb- 
ances, aurorae, and various irregular distributions ionospheric electrons. 
shows that these may all related simple way single causal agency, 
namely, certain convection system the outer portion the earth’s magneto- 
sphere. The source this convection taken viscous-like interaction 
between the magnetosphere and assumed solar wind, though other sources 
equivalent nature may also available. The model capable accounting 
for many aspects the phenomena concerned, including the morphology 
auroral forms and the occurrence ‘spiral’ patterns the loci maximum 
intensities several features. also bears directly the steady state the 
magnetosphere, and particular the production trapped particles the 
outer Van Allen belt. short, provides new basis which full under- 
standing these several phenomena may time built. 


INTRODUCTION 


Throughout the whole the earth’s exosphere, and through much the 
underlying ionospheric region, the geomagnetic field exerts strong control 
the motion ionized matter. This composite region geomagnetic 
dominance conveniently termed the magnetosphere, following Gold (1959). 

common think the greater part the magnetosphere static 
terms, relatively calm expanse quiescent plasma. But, Gold has 
emphasized, there real reason for adopting such view advance and 
indeed convective motions the ionization can established without diffi- 
culty. Convection the whole magnetosphere fact implicit variety 
circumstances (including, for example, those described the tidal dynamo 
theory) which driving forces are exerted the ionization relatively low 
levels. The complementary case convection generated extreme altitudes 
has not been paid much attention yet, but its consequences may impor- 
tant many phenomena. 

examine here one particular pattern convection due distant forces, 
whose natural occurrence are led expect. The driving mechanism that 
envisage consists viscous-like interaction between the outer magneto- 
spheric material and the interplanetary gas beyond. similar convection 
could established other processes, and one such proposed independently 
accompanying paper Fejer (1961). Regardless its manner genera- 
tion, however, the convective system discuss has consequences far- 
reaching import, and these that wish emphasize the present 


paper. 
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The convection system will relatively weak times low solar activity, 
but even then provides single comprehensive model explain the high- 
latitude patterns occurrence magnetic agitation, sporadic-E and spread-F 
ionization, radio absorption, and radio-star scintillation centers. With 
strengthened convection times high solar activity, the model can account 
for many features auroral excitation, morphology, and motions, and 
related magnetic variations. provides, new basis which theories 
detail all these subjects may time built. 


GENERAL DEVELOPMENT 


The Basic Model 

The outer boundary the magnetosphere all probability determined 
some way interaction between the geomagnetic field and the inter- 
planetary gas that lies beyond. have yet certain information the 
density and motion the interplanetary gas, however, our picture the 
outer magnetosphere must tentative best. For purposes the present 
discussion, shall adopt model which the interplanetary gas generally 
motion outward from the sun ‘solar wind’ (Parker 1958a), although 
shall assume that the strength this wind varies between quiet and dis- 
turbed conditions. shall ignore the motion the earth its orbit, although 
note this point that such motion could duplicate the essential feature 
(the viscous-like interaction) that wish abstract from the solar-wind 
model, even the model should prove untenable quiet conditions. 

must concerned briefly with the distribution geomagnetic field 
lines within the magnetosphere, especially the case the ‘high-latitude’ 
field lines that traverse the polar ionosphere. not known yet whether 
these field lines extend indefinitely into the interplanetary gas whether 
they form closed loops from one polar region the other, although the dis- 
tinction some respects fundamental. The model adopt here one 
which the high-latitude field lines indeed form closed loops. These loops 
not lie symmetrically about the geomagnetic axis, but instead are confined 
‘geomagnetic tail’ that formed the ‘lee’ the earth the viscous 
effect the assumed solar wind. This model has been described its simplest 
form Johnson (1960), while Piddington has discussed similar 
model for periods geomagnetic storms and Beard (1960) has studied some 
details its formation. Its topology implicit the model developed 
Chapman and Ferraro for the initial phase geomagnetic storms (Chapman 
and Ferraro 1931) and explicit Dungey’s model rotating magnetic stars 
(Dungey 1958). sketched here meridional section (Fig. and 
equatorial section (Fig. following schematically the form depicted 
Johnson. (The geomagnetic and geographic axes will assumed coincident 
and perpendicular the ecliptic plane for simplicity.) 

will noted that, this model, there bifurcation field lines 
the noon meridian: low-latitude lines are roughly dipole shape whereas 
high-latitude lines sweep round over the poles into the geomagnetic tail. 
Indeed, discussions the whole magnetic field convenient think 
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Fic. section the magnetosphere, taken the noon geomagnetic meridian. The 
solar wind (which directed away from the sun) moves from left right the diagram, 
causing high-latitude field lines swept into the lee the earth form geomagnetic 
‘tail’. The direction rotation the magnetosphere indicated. 


two classes field lines, the one lying low-latitude torus and having 
dipole-like shape, the other lying the high-latitude torus that constitutes 
the geomagnetic tail. Field lines that rise the day side the earth near 
the boundary that separates these tori will markedly deformed, and their 
allocation one other torus will means clear-cut, while the 
night side the distinction between the two tori essentially disappears. region 
confusion may defined the grossly distorted field lines that pass near 
the outer limits the geomagnetic field, and this region maps along the field 
lines into two roughly semicircular ‘zones confusion’ the earth’s surface, 
centered the noon meridian. These zones are depicted the accompanying 
figures. 

company with many others, and already noted, envisage solar 
wind that varies strength with solar activity. the wind increases, more 
and more field lines are likely carried into the geomagnetic tail, and the 
zone confusion would then proceed lower latitudes. This behavior 
reminiscent the behavior the auroral zone times high activity, and 
there very likely physical relationship between the two. Nevertheless, 
resist the tendency associate the zone confusion with the auroral zone 
itself, and instead place somewhat higher latitudes. For example, our 
subsequent illustrations will drawn the basis transition from low- 
latitude dipole field high-latitude magnetic tail latitude 70°, 
some above the conventional auroral zone. This transition corresponds 
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SOLAR WIND 


ZONE OF 
CONFUSION 


Fic. equatorial section the magnetosphere viewed from above the north pole, 


showing streamlines for the case which the magnetosphere simply rotating the manner 
depicted Fig. The streamlines are also equipotentials the electric polarization that 
induced the rotation. Neighboring equipotentials are separated 


cutoff the magnetic field geocentric distance earth radii 
the sunward side, using unperturbed dipole extrapolation, somewhat 
closer cutoff more self-consistent model. intended represent 
moderately disturbed conditions, produced solar wind perhaps 
proton-electron pairs per cm* moving speed 1000 km/sec. 


The Impressed Convection 

now assume that the interaction between the outer magnetosphere and 
the solar wind contains viscous-like component. this mean simply 
that some the momentum the solar wind transferred across the boundary 
the magnetosphere the ionization within. The nature this momentum 
transfer is, for present purposes, minor importance; its existence, the 
existence equivalent mechanism, crucial. 

ourselves tend think terms coupling initiated instabilities 
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the boundary the magnetosphere and carried depth form 
hydromagnetic eddy viscosity. The existence the requisite instabilities has 
been anticipated theoretically Dungey (1955; 1958, pp. 151-152) and 
Parker while Sonett (1960) has discussed observational evidence 
the subsequent coupling they produce. Sonett infers appreciable exchange 
momentum throughout the geocentric range earth radii during rela- 
tively quiet conditions. Our subsequent illustrations indicate turbulence over 
the range earth radii the low-latitude torus for moderately disturbed 
conditions. 

worth noting that other processes might produce the convective motions 
deduce here, but that our later discussion essentially independent 
the manner which the convection established. The alternative processes 
include, for example, the charge separation deduced Chapman and Ferraro 
(1933), the hydromagnetic interactions outlined Dungey (1961), and the 
asymmetric injection energetic charged particles now proposed Fejer 
(1961). fact, any mechanism that produces the pattern magnetic 
variations consequence electromotive forces generated above the 
ionospheric region must produce pattern magnetospheric circulation 
the form shall discussing here. 

the presence viscous interaction, which now assume for purposes 
discussion, the solar wind its passage past the earth will carry material 
the outer regions the magnetosphere away from the sun, towards the 
geomagnetic tail. This material, being linked the geomagnetic field, con- 
strained remain within the magnetosphere and builds excess 
concentration the region the tail. The build-up cannot persist steady 
state, and consequently return flow ionization must set the 
interior the magnetosphere. Circulatory loops will established, 
form depicted Fig. which the ionization convects continuously. the 
popular terminology hydromagnetics, the ionization will ‘frozen’ the 
geomagnetic field lines (see Dungey 1958, for example) and will carry them 
along the convective motion. The field lines can distorted only slightly 
the process, except perhaps the outer boundaries, and hence the whole 
each field line will carry out the convective motion imposed upon it. The 
ionization lower heights, being frozen onto the field lines, will similarly 
convect. Hence the pattern motion whose equatorial section shown 
Fig. will extended throughout the body the magnetosphere along the 
field lines. this gigantic convective system that forms the central feature 
our subsequent discussion. 

The foregoing (hydromagnetic) description may not satisfy all tastes, and 
shall rephrase accordingly more classical terms. begin again with 
the basic assumption that the ionization the outer magnetosphere set 
into motion the direction flow the solar wind, and that establishes 
certain steady-state velocity pattern (V). This motion the electrons and 
ions takes place the presence the geomagnetic induction (B), however, 
and can therefore proceed only electric polarization field (E) estab- 
lished the accumulation space charge offset the Lorentz force. (In 
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SOLAR WIND 


Fic. motion impressed the magnetosphere viscous-like interaction with the 
solar wind. This equatorial section the magnetosphere with the solar wind blowing 
from top bottom. The viscous-like interaction causes tubes force which lie near the 
surface the magnetosphere pulled around into the geomagnetic tail, indicated 
the large arrows. return flow (indicated the small arrows) takes place the interior 
the magnetosphere. 


the models Chapman and Ferraro (1933) and Fejer (1961) the process 
charge separation more fundamental, but would induce turn the 
motions now under discussion.) The field must derivable from 
potential (¢) the steady state, whence 


(Rationalized m.k.s.a. units are used throughout the formal development.) 
This relation limits the type velocity field that can established the 
magnetosphere, for implies that must everywhere lie the local equi- 
potential surface, and that that surface must oriented include the 
local vector. 

Consider now section cut through the magnetosphere surface orthog- 
onal the (magnetic) equatorial plane, for example. The equipotential 
surfaces will cut this surface some pattern closed curves which are them- 
selves equipotentials. Since the vectors lie the equipotential surfaces, 
their components (V,) orthogonal must lie along the equipotential curves 
and must trace out flow lines that are closed convective loops. The simplest 
and most readily justifiable pattern such loops consistent with the impressed 
motion the form presented Fig. which therefore recover the 
present reasoning. This pattern will examined more detail shortly. 
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Since equation (1) implies that the geomagnetic flux lines are equipotentials, 
also implies that the electric fields existing one region must mapped 
along the flux lines other regions. then implies that the magnetospheric 
ionization those regions also convecting. Indeed, not difficult 
show from equation (1) and relations that all the ionization that 
lies along tube magnetic flux one time must convect together, 
remain always along some common tube flux. This the classical analogue 
(part of) the ‘frozen field’ concept hydromagnetics, and recovers the 
earlier conclusion that the motion ionization depicted Fig. will 
mapped throughout the magnetosphere along magnetic flux lines. 

Certain factors tend inhibit the magnetospheric motions, 
warrant brief comment. the first place, energy would required given 
when change made from one steady state another. This energy 
would have transferred across the magnetospheric boundary, and 
presumably transferred when, for example, strengthening solar wind 
blows more and more field lines into the geomagnetic tail. But the steady 
state, the convective motions that are interest here contain far less energy 
than does the geomagnetic field. Accordingly, any field deformation that 
produced the convection itself should negligible (except perhaps the 
outermost levels, already noted, where the convective energy density 
might comparable with the geomagnetic energy density). hydromagnetic 
parlance, this conclusion implies that tube magnetic flux must vary its 
cross section convects, order maintain the strength the field 
its appropriate equilibrium value. The ‘frozen in’ ionization must similarly 
alter its cross section, and two-dimensional divergence will 
produced. (This latter statement may verified classical grounds the 
application equation (1) and relations once again.) This intro- 
duces strong tendency for the magnetospheric matter vary density 
and pressure convects, despite the possibility ionization ‘leaking out’ 
into the lower ionosphere, and such variations raise the question stability. 

The foregoing points have been stressed Gold (1959) extensive 
discussion the principles magnetospheric convection. the matter 
stability, notes that certain regions (defined shells magnetic field 
lines) might permit rapid adiabatic convection, while others, unless strongly 
forced, might permit only slow isothermal convection; still others might 
highly stable against any convective processes, and the Van Allen belts may 
particularly important this regard. suggests, however, that the 
outer magnetosphere likely such nature permit the more 
rapid convection, and that this convection might affect certain geophysical 
the very ones, fact, that will come under close study our 


features 


subsequent development. assume here that stability will not major 
problem the outer regions that are primary concern the present study. 

Atmospheric viscosity and ohmic losses become appreciable near the base 
the ionosphere and act inhibit the motion. Being purely resistive 
nature, they cannot prevent the convection but can only limit its speed. 
must assume that the observed velocities are determined the balancing 
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energy transfer the magnetospheric boundary with energy loss these 
dissipative processes, once any inherent stability has been overcome. This 
assumption will bear further analysis another time, but two preliminary 
checks can made quite readily. The losses are mainly ohmic, and 
amount something like watts/m for typical auroral currents. This loss 
extends throughout height range the order 104 over area the order 
for total dissipation rate the order watts. The local rate 
energy dissipation down order magnitude from the rate introduced 
normal photo-ionization processes, which comforting, while the total 
rate down order magnitude more from the rate energy capture 
implied Sonett’s (1960) calculations. The uncertainties inherent both 
the present estimates and Sonett’s leave open the possibility that even 
closer agreement could obtained imply that the greater part the energy 
incident from the solar wind fact dissipated ionospheric currents with 
magnetospheric convection acting intermediary. Such conclusion 
would overcome the difficulty encountered Sonett the search for sink 
for the energy influx deduced. 

this point, might add some credence our views viscous-like 
interaction the following comments. shall later argue that the motion 
auroral irregularities corresponds the motion associated geomagnetic 
field lines (using the hydromagnetic terminology). The observed speeds are 
normally somewhat less than km/sec, and these would scaled 
something like km/sec for the speed the corresponding field lines their 
greatest height above the earth. This speed still substantially below the 
speed associated with the solar wind disturbed times, must the 
operating mechanism indeed viscous-like. 

Turning now the details the convection pattern, note that these 
will depend large measure the stability and the dissipative processes 
discussed above. For example, the return flow within the body the magneto- 
sphere would tend channelled through regions low rather than high 
stability. Energy dissipation ionospheric heights results from current flow, 
and the currents must driven forces the magnetosphere. These forces 
would include not only the basic interaction with the solar wind, but secondary 
mechanisms such those provided ion and electron pressure gradients. 
The net steady-state flow would determined the interplay these 
several features, and would governed formidable array relations 
which equation (1) only single member and approximate one that. 

shall make attempt here any analytical approach this difficult 
problem. Instead, shall simply sketch refined version Fig. form 
which seems reasonable us, which based large part observational 
evidence (assuming our interpretation the observational material 
true), and which compatible with the physical principles the problem. 
The proposed pattern presented Fig. Certain its features warrant 
discussion. 

Our subsequent comparison with observational data leads place the 
innermost part the return flow geomagnetic shell whose position (as 
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SOLAR WIND 


Fic. The proposed pattern streamlines (or alternatively the equipotentials the 
electric field, separations kilovolts) the equatorial plane the magnetosphere. The 
main circulation does not penetrate closer than about earth radii; dotted lines within this 
region indicate possible inner convection system. 


measured the base the field lines) somewhat below the conventional 
auroral zone, latitude about 62°. This corresponds return flow 
bounded its inner side geocentric distance earth radii the 
equatorial plane. These values are probably representative moderately 
disturbed conditions and some suitable mean between lower numbers 
times strong disturbance and higher numbers during relatively quiet 
periods. 

would associate the inner boundary region with some form stability, 
and consider most likely that the energetic particles the outer Van 
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Allen belt provide the stabilizing agency. The suggestion that this belt 
found lower latitudes the night-side the earth (Rees and Reid 1959) 
would imply turn that the circulatory system could penetrate greater 
depths and lower latitudes night. similar conclusion might drawn 
for other potential stabilizing processes, and would help explain the 
morphology aurorae. is, however, refinement that shall not attempt 
incorporate our illustrations. 

may also note passing that the circulatory system will some extent 
self-stabilizing, for results the energization particles during the 
inward convection will discussed. The energization tends increase the 
longitudinal drift motions that are introduced the inhomogeneity the 
geomagnetic field, and these motions then tend carry the convecting particles 
around longitude out the zone inward convection. 

addition the main convective patterns already discussed and depicted 
solid contours Fig. note the possible existence inner circulatory 
system represented broken contours. This inner circulation complicates 
the suggested pattern manner analogous certain simply hydrodynamic 
systems, and might set some viscous (hydromagnetic) interaction 
the inner boundary the outer circulatory system some other 
secondary process such will described subsequently. postulated here 
primarily for the explanation the low-latitude system currents, 
although may have other useful implications well. will discussed 
only briefly the subsequent development. 


Superimposed Effects Rotation 

The discussion this point has treated the convection departure from 
static condition. fact, course, the earth and its ionosphere are rotation, 
and the magnetospheric ionization will also rotate (as consequence hydro- 
magnetic coupling) provided drag effects imposed the interplanetary gas 
are not too great. The rotation the high-latitude torus takes place about its 
own axis, rather than the earth’s, Johnson (1960) has described, and 
freed from the inhibiting effects that would have arisen (Hines 1959) had the 
high-latitude field lines been taken extend into the interplanetary gas. 
the present model, then, seems reasonable assume that rotation 
proceeds without significant impediment except perhaps the very boundaries 
the magnetosphere where the viscous interaction imposed. 

with the convective motion, the rotational motion can proceed only 
the presence electrostatic field. The equipotential surfaces the field 
will again lie along the magnetic flux lines and the rotational flow lines. Their 
equatorial sections have been depicted already Fig. 

There must some interaction between the rotational motion and the 
convective flow, but its evaluation lies beyond the scope the present study. 
shall content with derivation the net flow produced simply 
superimposing the two independent velocity fields. The superpositioning can 
achieved most simply treating the contours Figs. and equi- 
potentials, assigning values potential individual curves each set 
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contours, adding the two values potential number points throughout 
the field, and drawing new contours constant total potential. The assign- 
ment values the original contours straightforward the case the 
rotational system, least for the undistorted portions the geomagnetic 
field, since the angular velocity and the field are known. The corresponding 
assignment for the counter-rotating geomagnetic tail and for the convective 
system rather arbitrary the absence detailed analysis, although the 
direction increasing potential predetermined the direction required. 
Despite this arbitrariness, the general trend the resultant flow pattern 
quite unambiguous; depicted for one specific case Fig. The principal 
feature noted that the ionization convecting inward the nighttime 


SOLAR WIND 


Fic. composite picture obtained superimposing the equipotentials Figs. (2) 
and (3), showing the type motion expected the presence rotation and viscous- 
like interaction between the solar wind and the surface the magnetosphere. 
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side the earth carried the rotation around and through the morning 
side. Other features will pointed out later, they become pertinent the 
discussion. 


The Motion Energetic Magnetospheric Ionization 

The convective motion described far has taken account the gyrations 
individual particles about the geomagnetic field lines their random 
motions along those lines. Such motions lead additional drift components, 
now well known connection with so-called trapped radiation. The 
additional drift will not important for particles sufficiently low energy, 
but becomes quite comparable with the convective motions for energies 
the order few kiloelectron volts. Particles greater energy still will 
move predominantly with their natural drift velocities, and will subject 
only small perturbations from their natural orbits consequence the 


SOLAR WIND 


CONVECTION 
ENERGETIC ENERGETIC 
ELECTRONS PROTONS 


Fic. effect the proposed circulation pattern energetic trapped particles. 
thin shell such particles split into separate electron and proton shells, such that the 
electrons drift lower-latitude field lines during the morning hours and the protons during 
the afternoon hours. 
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magnetospheric convection (or, more properly, consequence the electric 
fields that accompany that convection). Since proton and electron drifts 
proceed opposite senses around the earth longitude, while the convective 
effects act the same sense any one point, given shell (mono-) energetic 
protons and electrons will split into two deformed shells, one protons 
and one electrons, depicted Fig. 

clear that this effect will tend produce charge separation the 
energetic particles, even when integrated over range energies and through 
continuous distribution shells. Fejer (1961) points out that this separation 
will modify the convective system, and that will such way 
produce the inner loops tentatively drew Fig. That figure should 
now considered composite, including the primary convective system 
and any superimposed secondary effects the charge separation. The latter 
tend only transitory, however, since steady state the separation 
energetic charged particles would largely offset the accumulation 
low-energy space charge; some residual motions, having complicated form, 
would remain. 


The Induced Ionospheric Motions 

The convection systems that have far been depicted equatorial 
sections take somewhat different appearance when extrapolated down 
the geomagnetic field lines. The equipotential surfaces have sections iono- 
spheric heights that are best illustrated polar view the earth, 
Fig. for the convection alone and Fig. with rotation superimposed. The 
outermost regions the magnetosphere map into the ‘zone confusion’ 
previously introduced, and this zone shown explicitly. The flow pattern 
and near the zone difficult deduce because the severe distortion the 
pertinent field lines. have contented ourselves the present representation 
with mapping the flow pattern the more certain regions distance 
from this zone, and with implied continuity flow across it. Figures 
and show the equipotentials accordingly. Residual deviations 
pattern may have some bearing geophysical observations, however, and 
should kept mind. 

the ionospheric region, equation (1) remains valid high degree 
accuracy, and continues imply convection ionization along the equi- 
potential surfaces. The same true for electrons the region, down 
heights (say) 100 km. the case ions, however, collisions become 
important and must taken into account. There results deviation the 
ion motion from the simple equipotential contours. 

The effect perhaps best seen from the equation 


which better approximation than (1) the equation motion ions. 
(Here the charge ion, its mass, and its velocity; the 
velocity the neutral gas, and ‘frictional frequency’ related the 
mass neutral molecule. Inertial and certain other terms are still missing 
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from the equation given.) may think first the case when and 
then note that (1) likely good approximation (2) only ifeB/M,>K,, 
that is, the ion gyro-frequency greatly exceeds the ion frictional frequency. 
The two become comparable the upper part the region, however, and 
the lower part this region the gross inequality fact reversed. the 
latter circumstances, approximates and considerably reduced 
magnitude from the value attained ions higher levels and 
electrons even the region. The diminished velocity relative electrons 
implies that the electrons will become the main current-carrying constituent, 
and that the current will then flow roughly along the equipotentials but 
the sense opposite the motions indicated Fig. (Figure appropriate 
this connection rather than Fig. because the rotational component 
the collisional term equation (2) offset equal component 
there, and (2) may replaced once again (1) insofar the rotation 

These conclusions require more exact analysis terms the anisotropic 
conductivity tensor that applies E-region heights, and with the inclusion 
resultant polarization fields. The net effect would 
effect, and its nature would depend height-integrated conductivities. The 
latter have been deduced various authors (for example, Baker and Martyn 
1953) for application atmospheric-dynamo theories, and has been found 
that the Hall conductivity predominates except equatorial latitudes. This 
conclusion leads, effect, the description just given which the electrons 
are the dominant carriers current. 

The development polarization fields certain circumstances will lead 
flow current regions not directly affected the convection process, 
and this secondary current can turn generate secondary convection systems. 
This complication provides another possible source the internal circulation 
depicted Fig. and the associated low-latitude current system. 

Unlike atmospheric-dynamo theories, which the motion the neutral 
atmosphere prescribed, the present theory involves the neutral gas with 
convective motions that are impressed ions. These may signi- 
ficance under certain conditions. They are most simply studied application 
the approximate equation motion 
(3) 
which the driving force taken solely the reaction the neutral gas 
the collisional term (2); and are the number densities neutral 
particles and ions respectively. The implications this equation are turn 
most simply studied under the assumption that and vary exp 
and that d/dt then 


where 
(5) tk = Na 


evident from this that the convective motion the neutral gas likely 
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negligible only for variations time scale much shorter than 
which the order some tens minutes the region and some hours 
the region. 

For variations with time scale longer than the other hand, 
evident that the neutral molecules can carried along the ions their 
motion, provided other forces interfere. The residual difference velocities, 
for wrx given 


(6) —twtrV, 
from (4), and this result combines with (2) yield 


the approximate equation motion for the ions. now see that the ions 
will follow the convective motion (1) provided eB/M, which 
less stringent condition than the earlier requirement that eB/M, Indeed, 
evident that the ions may convect closely with the electrons motions 
whose time scale least long 


(8) 


even heights where the ion collision frequency exceeds the ion gyro- 
frequency. The quantity varies throughout the region and varies from 
day night; the order few hours 125 during the day, and 
increases from this value night and lower levels unless abnormal 
sources ionization are introduced. are led expect the duration any 
intense current flow limited this effect. Residual currents will flow 
for times exceeding but they will less intense and their pattern may 
differ appreciably from that the equipotentials. These aspects will 
discussed again relation the observational data. 


The Distribution Turbulence 

has been noted already, expected and observed that the outer 
portions the magnetosphere are turbulent. Our present model indicates that 
the material these regions will convected the interior, and therefore 
expect corresponding convection the turbulence. The turbulence energy 
will course decay the convection proceeds, but this may more than 
offset the effects the compression that accompanies transport the 
interior. Such compression would result increased absolute fluctuations 
ionization density and magnetic field, the fractional fluctuations remained 
essentially unchanged, and increased gradients due the diminishing linear 
dimensions. The way then open for the acceleration particles high 
energies (Parker 1961), and for the production severe hydromagnetic 
disturbances. clear that variety geophysical phenomena might follow, 
and the pattern convected turbulence therefore some consequence. 

This pattern depicted Fig. which retains the major convective system 
Fig. and introduces the following regions for convenience discussion: 

The outer magnetospheric regions which the turbulence originally 
generated. Here the turbulence exists independently the circulatory motions, 
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SOLAR WIND 


distribution turbulence the equatorial plane the magnetosphere. The 
hatched area denotes the region which there intimate contact with the interplanetary 
gas and into which momentum transferred directly the solar wind. The various zones are 
the text. 


and may expected decrease intensity from the outer boundary 
inwards. 

The regions through which the outermost material convected. 
This material likely the most turbulent initially, and subject 
strong compression. may then expected major concern. 

The region between and the morning side. The material here 
convected from the interior and not severely compressed. should 
therefore turbulent, but not the same importance 

and Regions containing turbulence convected from the interior 
the evening side. loops round the earth the morning side, and the 
evening side, although will not exist unless the convection sufficiently 
strong relation the rotation. The material not turbulent 
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that initially, but subject stronger compression; therefore 
ranks with its possible geophysical consequences. The material 
even less turbulent initially and subject slightly less severe compression 
than that IV, but the compression still strong and its turbulence may 
therefore significance those occasions when the region itself exists. 

VI. region that essentially devoid turbulence, since the material 
contained never convects through the convection sufficiently 
weak, that does not exist, extends from the inner portion the 
magnetosphere the evening side. 

VII. The inner region into which the circulation does not penetrate. 
Turbulence generated the surface the magnetosphere unimportant 
this region, which corresponds geomagnetic latitudes less than about 62° 
the surface the earth. 

The regions depicted Fig. map along the field lines into specific zones 
ionospheric heights and ground level. These are illustrated Fig. 10, where 
the numbering system Fig. retained. The ‘zone confusion’ previously 
introduced lies centrally along zone this figure. 

view widespread interest ‘spiral’ patterns high-latitude geo- 
physical phenomena, worth noting even this stage the appearance 


NOON 


MIDNIGHT 
The distribution turbulence the magnetosphere observed ionospheric 


levels. This pattern obtained mapping the distribution Fig. down along lines 
force the geomagnetic field. 
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spiral-like forms zones II, IV, and This feature will receive more dis- 
cussion directly, but quite clearly has relation the ‘Stérmer spirals’ 
that are often invoked. The latter can have little direct bearing the 
observed phenomena any event (Agy 1960). 


GEOPHYSICAL APPLICATIONS 


Features Associated with Magnetospheric Turbulence 

Whenever sufficiently intense turbulence occurs the magnetosphere, its 
effects should observable means magnetic and ionospheric measure- 
ments ground level. particular, might expect such diverse phenomena 
geomagnetic agitation, auroral activity (both visual and radio), auroral 
absorption, high-latitude types sporadic-E and spread-F ionization, and 
high-latitude radio-star scintillations associated part with magneto- 
spheric turbulence. Indeed, the field-aligned irregularities that characterize 
several these phenomena can perhaps best explained consequence 
just such turbulence, mapped along the geomagnetic field lines particle 
precipitation and hydromagnetic propagation. The details this mapping 
process postpone for separate investigation and turn instead compari- 
son the theoretical distributions and motions with readily available obser- 
vational data. 

some cases, the distributional information has been summarized the 
form geographic loci maximum intensity the feature under considera- 
tion. While such summaries run the risk being misleading, particularly 
when more than one maximum found, they are nevertheless convenient and 
adequate for our immediate purpose. Summary data this type are displayed 
Fig. for geomagnetic agitation (Burdo 1957; see Hope 1961), intense 
high-latitude ionization (Hagg al. 1959, December data), 
ionization (L. Petrie, private communication), visual aurorae (Malville 
1959), ‘auroral-type’ radio scattering (Forsyth 1960), and auroral 
absorption (L. Montbriand, private communication). 

readily apparent from this display that major enhancement 
disturbance phenomena found single broad swath that sweeps 
from higher lower latitudes over range longitudes the night side 
the earth. This duplicates form our combination zones and 
their compression phase, and our earlier arbitrary (but justifiable) choice 
boundaries for the magnetospheric convection was such ensure duplication 
the latitude range well. are led conclude that the dominant dis- 
tribution high-latitude disturbance phenomena can explained adequately 
consequence convected turbulence, following the pattern given our 


model. 
The two secondary spirals provided the observational data the morning 


and afternoon sectors are not readily explained. The characteristics 
magnetic agitation these spirals are sufficiently different from those the 
nighttime loop that some authors have sought separate mechanism for their 
source (for example, Hope 1961). Comparison with Fig. then leads 
suggest that they are associated with zone and hence with the primary 
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NOON 


MIDNIGHT 


Fic. The positions diurnal maxima various disturbance phenomena geomagnetic 
co-ordinates. The loci and refer geomagnetic agitation and visual aurora; 


unconvected turbulence the outer magnetosphere, perhaps combined 
some extent with the recovery phases zones II, and This view cannot 
strongly supported, however, prior fuller examination both the 
theory and the observations. 

Scintillations radio stars, due apparently ionization irregularities high 
the region, have not been analyzed the fashion Fig. 11. But observa- 
tions made subauroral latitudes, related ionization low auroral 
latitudes, reveal peak scintillation index near geomagnetic midnight 
(Hewish 1952; Dagg 1957). The index rises rapidly the late evening hours, 
and falls gradually during the morning. This behavior appears consistent 
with the dominant role the combined zone, The earlier onset 
scintillations times high activity (Dagg 1957) would follow naturally 
upon the expansion the circulatory loops Fig. lower latitudes, and 
upon the increased role that might played zone 

Certain evidence supports the suggestion that there region low dis- 
turbance, namely the non-turbulent zone VI, situated the evening side 
the auroral zone. For example, representative mean diurnal variation 
ionospheric absorption has been published for College, Alaska, geomagnetic 
latitude 65° (Little and Leinbach 1958). While does not depict any clear 
localized maxima, does reveal very marked minimum absorption some 
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hours before geomagnetic midnight. The diurnal variation sporadic-E 
Narssarssuak (geomagnetic latitude 71°) shows broad and pronounced 
minimum occurrence this time (Olesen and Wright 1961). The non- 
turbulent zone would course also associated with the quiet auroral arcs 
that characterize the pre-midnight period contrast the much more active 
forms the early morning hours (for example, Davis 1960). 

Our picture turbulence leads expect that any given 
magnetospheric irregularity would convect more less with the pattern flow 
depicted Fig. The corresponding ionospheric irregularities would then 
follow the pattern shown Fig. or, when measured relative the rotating 
earth, that shown Fig. The available evidence supports this conclusion. 

The evidence derives part from measurements radio-star scintillations. 
The ionization irregularities, which lie low auroral latitudes the cases 
interest, move predominantly westward before geomagnetic midnight and 
eastward after (Hewish 1952; Briggs 1960). The tendency for reversal 
midnight becomes more pronounced the magnetic activity increases, which 
suggests that the convective systems have been discussing may com- 
peting with some other mechanism (for example, the tidal dynamo) the 
relatively low latitudes concerned. Other evidence this effect can adduced, 
but real importance here. 

Much more complete observations are available from auroral studies. 
Measurements made low auroral latitudes reveal again the tendency for 
westward motion before midnight and eastward motion after, both 
visible and radio aurorae (for example, Meinel and Schulte 1953; Kim and 
Currie 1958; Bullough and Kaiser 1955; Lyon and Kavadas 1958). has also 
been established explicitly (Kim and Currie 1958) that systematic westward 
motion introduced the eastward rotation the earth, and this conclusion 
would difficult explain unless eastward rotation were associated with 
the generating mechanism. Such rotation is, course, inherent the model 
have developed, and the observational evidence compatible with the 
implications Fig. 

The pattern motion auroral irregularities has been established more 
fully Davis (1960), analysis IGY data. His work includes observa- 
tions made high auroral latitudes, where the east-west motions are found 
reverse just the motions depicted Fig. reverse. Indeed, the over-all 
pattern motion presented Davis essentially duplicated the pertinent 
portions Fig. Moreover, the alignment auroral forms the 
direction geomagnetic latitudes greater than about 80° would suggest that 
the flow lines for this region shown Figs. and are essentially correct, if, 
lower latitudes, the auroral forms tend lie along flow lines. (The 
explanation this tendency lies beyond the scope the present paper, 
does the explanation auroral forms themselves. suspect that both follow 
from form instability the associated currents, and this possibility 
now being examined more detail.) Davis’ results have been recently con- 
firmed Denholm (1961) and Stoffregen (1961). 

The description the motions aurora and related phenomena can 
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completed allowing for additional motions. have suggested, 
for example, that the outer radiation belt prevents the circulation from 
penetrating deeper than about earth radii the equatorial plane the 
magnetosphere. Due longitudinal variations the strength the geo- 
magnetic field, the particles forming the belt tend move lower latitudes 
the night side the earth (Rees and Reid 1959). This effect can expected 
result motion the low-latitude boundary auroral 
activity whole, the motion being towards the equator the evening and 
away from the equator the morning. Such motions have been described 
practice Bond (1960) for example. This movement, which may amount 
several hundred kilometers, best observed from stations outside the auroral 
zone; should not confused with the motion irregularities 
within auroral form, which not necessarily follow the same pattern. 

Finally, might note obvious fact: the increase the speed 
scintillation and auroral irregularities that accompanies increased magnetic 
activity (for example, Dagg 1957) immediately explained the present 
model, due simply more rapid convection induced strengthened 
momentum interchange the magnetospheric boundary. The latter could 
produced increase either the strength the irregularity the 
solar wind. 

summarizing all these points, may said that convected turbulence 
clearly capable accounting for the distribution and motion wide 
variety high-latitude geophysical features. believe that such explana- 
tion likely prove correct. 


Features Associated with the Large-Scale Convection and Compression 

addition their apparent enhancement turbulent features, magneto- 
spheric convection and compression are capable more direct effects 
larger scale. the case convection, these can come about the displace- 
ment energetic particles already trapped the geomagnetic field (as 
described and the transport energetic solar particles from the 
boundary the magnetosphere its interior. The compression, which 
fact inherent aspect convection inhomogeneous field, corresponds 
energization both types particle. 

The energization process derives from motion ions and electrons across 
the equipotential surfaces and therefore implies violation equation 
(1). The cause this violation found the drift motion that ionization 
naturally attains the inhomogeneous geomagnetic field, and this motion 
proportional the energy the ionization amongst other things. Whether 
the drift speed large small comparison with the convective speed then 
depends the energy: low-energy ionization tends follow the flow lines 
the convective system, while high-energy ionization tends follow its 
natural drift trajectories. transitional energy (W,) may envisaged 
which the two components motion are comparable magnitude. Its value 
varies with the strength the convective system and with geocentric distance, 
but lies typically the range kev. 
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For the lower-energy ionization particular, often useful consider 
the energization process one adiabatic compression (Gold 1959). 
losses through precipitation are ignored, ionization convected from geo- 
through the diminution its cross section normal the (dipole) 
geomagnetic field, and further factor the order due the 
diminution the length field line above the precipitation level. Pre- 
cipitation into the lower ionosphere reduce these factors, but they provide 
preliminary basis for discussion. raised the power y—1 (where the 
adiabatic index, specific heat ratio thermodynamics) they yield the 
appropriate factors for the energization process. 

interparticle collisions are important, the effective 5/3 and the energy 
factor becomes collisions are unimportant, for the trans- 
verse compression and for the longitudinal compression (Spitzer 1956), 
the energy transverse thermal motions will altered the factor 
and that the longitudinal motions The convected low-energy 
ionization evidently alters its geocentric distance factor 
the course its circulation, its energy can increased order 
magnitude more. 

The higher-energy ionization cannot attain such large proportionate 
enhancement its energy, for the maximum ratio can attain 
limited its natural drift velocity discussed immediately above. this 
case, possible put limit the absolute amount energy that 
given particle can acquire, the collision-free case, since the energy derives 
essentially from the electrostatic field This maximum energy will 
the order eVBL, where representative length scale and the other 
parameters have been defined already. auroral heights 
whence W,, electron volts, where typically m/s and 
may taken roughly, and kev. This value might 
easily increased decreased order magnitude, depending the 
although the significance this conclusion has yet assessed. 

The particles the outer magnetosphere may expected exhibit 
broad range energies, from fraction electron volt solar-wind 
proton values several kev. The solar-wind protons are expected share 
their energy with electrons relatively short time Gold, private com- 
munication), thus providing supply electrons having energies the order 
kev. The convection will act all particles carry them inwards, 
although the more energetic will course break free they tend towards 
the energy The inwardly convecting stream will then act continuous 
source particles about this energy, throughout the night half the 
magnetosphere. Electrons this energy are capable penetrating auroral 
heights, and they may therefore looked potential source aurorae 
and abnormal ionization. 

The precipitation energetic particles fact depends upon reduction 
their mirror heights, and this takes place the particles are carried inwards 
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the night side the earth. The strength the magnetic induction (By) 
the mirror point particle whose velocity components the equator are 
and perpendicular and parallel the field respectively, given 


where the strength the magnetic induction the equator the field 
line interest. collisions are negligible, the ratio remains constant 
during compression. The resultant increased proportion the 
increasing energy the particle and the same time the field line which 
this must attained moving lower latitudes. Both effects lead 
depression the mirror height, which total may quite considerable. 

The foregoing effects are course more significant for the softer component 
the particle distribution (particle energies <20 kev, say) because the 
limit W,, placed the acquisition energy, and look this component 
the source most auroral forms. The harder component will not un- 
affected, however, and expect that important certain stable emissions. 
particular, suggest that the very stable band emissions that 
accompanies aurorae due the precipitation trapped protons 
moderately high energies (~100 kev, say) subject the effects just described. 
The emissions are found latitudes below the lowest-latitude 
auroral forms the pre-midnight period, and then move greater latitudes 
where they normally become lost the general auroral emission (L. Mont- 
briand, private communication). This crossover the sense expected 
from mapping the proton loop Fig. down auroral heights. 

perhaps worth noting that the convection and compression discussed 
here tend favor the precipitation protons the late evening and electrons 
the early morning. This consequence the fact that energetic protons 
drift westward, and energetic electrons eastward, from the region inward 
convection centered midnight (cf. Fig. 6). There observational evidence 
just such bias auroral emissions (Rees al. 1961). 

The precipitation both convected and trapped particles would course 
enhanced the action turbulence, with implications that have already 
been examined. But the mechanisms discussed the present section are not 
dependent turbulence, and can proceed even the non-turbulent region 
Fig. They can then responsible for the quiet auroral forms and 
glow that occur zone Fig. (and elsewhere), and believe them 
major contributing factor least. 

The fact that all our mechanisms—turbulence, convection, and compression 
tend favor the production aurorae night possibly too obvious 
warrant comment, although achievement not often attained auroral 
theories the past. 

not propose develop the fine points our model this stage, 
although easy see that the theory allows considerable scope for the 
understanding more detailed features than have yet discussed. would 
emphasize that our mechanisms auroral production are means 
exclusive, and would expect final description auroral phenomena 
include many aspects theories other than (but compatible with) our own. 
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Magnetic Disturbances 

The large-scale magnetic disturbances that contribute geomagnetic 
storm variations are often thought four parts: the initial-phase world- 
wide increase, the main-phase world-wide decrease, the diurnally varying 
component, and short-lived ‘bays’. Some writers consider the variation 
represent simply statistical averaging the bay disturbances, and its 
separate identity certainly has not been clearly established. 

There widespread agreement that the initial phase storm direct 
consequence solar particles approaching the geomagnetic field, the manner 
studied first Chapman and Ferraro (1931) and more recently recast 
hydromagnetic terms (for example, Dessler and Parker 1959). There also 
agreement that the main phase could caused the enhanced circulation 
trapped energetic particles (Singer 1957; Dessler and Parker 1959; Akasofu 
1960), the form ‘diamagnetic ring current’, although the means 
providing enhancement have been means clear. Piddington 
particular has challenged the processes far advanced for this purpose. 

direct application the development the preceding section indicates 
new process, and very powerful one that. The compression magneto- 
spheric ionization, and its consequent energization, would lead enhanced 
circulation the type required. This enhancement would offset some 
extent the decompression and de-energization that occurs the day side, 
and the loss magnetospheric ionization through enhanced precipitation; 
whether the balance would favorable not question that must 
bypassed for the present. But, addition the normal magnetospheric 
ionization, material captured from the solar wind will caught the 
convection previously described, and its contribution must the sense 
required for the main phase. therefore suggest that the convective system 
provides the mechanism introducing solar particles the interior the 
magnetosphere, and that supplying thereby the storm-time ring current. 
This mechanism avoids the difficulties discussed Piddington, whether these 
are real not. 

this picture, would interpret the delayed onset the main phase 
(by few hours, following the start the initial phase) the time taken 
before substantial solar material convected the interior and energized. 
The subsequent decay the main phase could follow from the capture process 
discussed Dessler and Parker (1959), which independent the initial 
generation process. 

Turning now the high-latitude system and bay disturbances, 
recall the current systems that would established convection 
discussed Section first approximation, these currents will 
flow along contours constant potential the sense opposite the ‘frozen’ 
convective motions electrons. Chapman’s analysis (1956) the system 
provides pattern current flow which essentially equivalent the 
convective pattern Fig. and with the direction flow reversed required. 

During storms greater intensity than those analyzed Chapman, 
appears that the diurnally varying current systems are shifted somewhat 
earlier hours and intensified auroral latitudes, strongly after midnight and 
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weakly after noon (Silsbee and Vestine 1942). These intensified currents are 
more often associated with ‘bays’ than with the system, but the semantical 
schism perhaps best avoided, when possible, until the phenomenon itself 
understood. 

advance two suggestions explanation the change, although 
recognize that the problem extremely involved detail and simple 
solution likely complete itself. the first place, note that the 
flow strong currents dependent the neutral atmosphere remaining 
free from the convective motion the ionization, discussed and 
that this freedom may lost after period few hours the convective 
forces are sufficiently strong. the neutral atmosphere the region rotates 
through the midnight position auroral latitudes, enters new circulatory 
loop which its motion not adjusted, and strong currents can flow. 
dawn, however, can have been accelerated take part the new convection, 
and the current strengths fall. The transitional sequence repeated the 
hours immediately following noon, when relatively strong currents flow once 
again. 

The second process applies particular the post-midnight enhancement 
and may account for its relatively great strength. involves simply 
enhanced conductivity due increased ionization, resulting from the combined 
effects the turbulent zone and the compression and convection 
discussed previously. 

both cases, and indeed the case the quiet system, divergence 
ionization flow must anticipated. This will result secondary system 
electric fields, and secondary current flow, such that the system 
Fig. has superimposed dipole-like current systems for which the direction 
the dipole that the undistorted local current. The nature these 
secondary currents has been discussed Martyn (1951), who has shown 
that the low-latitude system likely produced result. While 
not agree with Martyn’s method establishing the initial current flow, 
will discussed, consider the secondary electric fields important 
and capable producing the low-latitude current system. The secondary 
fields would turn modify the pattern convection have been discussing 
far, and particular they could establish convection the inner magneto- 
sphere the type indicated Fig. 

the basis this theory the current system should regarded 
being present long there solar wind and simply becoming greatly 
enhanced during geomagnetic storms. number effects may 
cause the pattern become (partly) apparent even during relatively quiet 
periods, such the enhancement conductivity result sudden iono- 


spheric disturbances. 

aside, may noted that Gold (1959) has suggested that low- 
latitude convection established normally, consequence strong 
equatorial heating. This would expected produce rising motion 
daytime, outward through the magnetosphere, and was suggested 
source low-latitude ionization irregularities. Without examining the merits 
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this view, would simply remark that such convection would opposed 
and inhibited the low-latitude convection just discussed, and that the 
occurrence low-latitude spread-F ionization indeed times 
geomagnetic disturbance (Lyon 1960). 

Returning now auroral latitudes, may note that the intensity and 
sense bay disturbances are related the speed and direction motion 
auroral irregularities (both visual and radio; for example, Heppner 1954, 
Harang and 1960). particular they exhibit positive negative 
changes the horizontal geomagnetic field according the auroral motion 
westward eastward. This association represents correlation detail 
the over-all statistical patterns already discussed, and confirmed 
even more detail the work Kaiser (1958) extremely short-lived bays. 
exactly what would anticipated the simple basis have adopted, 
which convection E-region heights largely confined the electrons 
alone. the case the short-lived bays, would probably have associate 
the enhanced speed motion with turbulence, but the same result would 
follow nevertheless. 

with the geophysical features previously examined, will not attempt 
develop our model detail present nor incorporate many attractive 
aspects other studies. Again evident, however, that the model does 
include properties whose further examination and amplification may help 
elucidate the details observed. 


The Equilibrium State the Magnetosphere 

have suggested that, even quiet times, circulation the type 
have discussed will present. so, should significance the equi- 
librium state the magnetosphere. Two examples will serve illustrate the 
role may play. 

First, consider the ionization thermal energies that constitutes 
far the largest part the material the magnetosphere. has been thought 
the past that the distribution this ionization above the region 
governed the hydrostatic law, suitably modified take account ambi- 
polar and rotational effects (for example, Johnson and Fish 1960; Opik and 
Singer 1960). However, the circulatory motion likely quite comparable 
magnitude with the thermal motion ions, and therefore likely affect 
the distribution appreciably. 

extreme case, which ion flow into and out the region was 
prevented, the ion content tube geomagnetic flux would remain un- 
altered convected, and steady equilibrium could established only 
that content were proportional the flux contained the tube. Near the 
equatorial plane this would result density roughly proportional 
where geocentric distance; other latitudes the distribution would 
more complicated functional form. This differs substantially from the 
hydrostatic result, and therefore incompatible with it. would expect the 
actual distribution determined appropriate combination the 
pertinent hydrostatic and convective forces. 
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turn now the more energetic ionization that constitutes the radiation 
belts, and consider particular the ‘soft’ component the outer belt 
which particle energies are the order kev. have already noted that 
energies this order can produced the convection times auroral 
activity, and even substantially weaker circulation could consequence. 
the circulation proceeded slowly and quietly, these particles would de- 
energized the course the subsequent outward convection, and they would 
make net contribution the trapped component. the presence tur- 
bulence, however, the pertinent adiabatic invariants break down, and 
substantial portion the energized particles can further energized and 
trapped shown Parker (1961). 

Parker stresses the point that turbulence the outer boundary the 
magnetosphere not sufficient, for more energized particles would diffused 
outward and lost than would introduced and retained the process. the 
other hand, strong turbulence the interior, surrounded relatively 
quiet region, would ideal for the production trapped particles. com- 
pressed turbulence effective our previous discussion implied, seems 
provide precisely the type situation required for the effective operation 
Parker’s process. Moreover, provides optimum conditions close the 
inner boundary the main convective system, geocentric distance 
4-5 earth radii which essentially the position the maximum flux the 
soft component the outer belt. are led suggest that this component 
due particles solar origin, convected the interior and there trapped 
the effects turbulence. envisage strong injection during severely 
disturbed conditions, but also lower level injection quiet times such 
will maintain the equilibrium flux density against losses the atmosphere. 

The variations that are found times magnetic storms also receive 
fairly natural explanation the present basis. The flux detected particles 
first shows marked decrease and recovers after some hours levels sub- 
stantially above normal (Arnoldy 1960). The preliminary decrease has 
been associated some with the general compression the geomagnetic 
field during the initial phase the storm, but should now note that the 
precipitation produced enhanced convection, and associated increase 
would provide alternative mechanism, and one which might well 
dominant. (If were, would expect the spectrum particle energies 
harden this time, since the convection process acts more strongly the 
soft component; but this conclusion might also follow from the suggested 
general compression.) this regard, important note that the enhance- 
ment the convection should proceed almost simultaneously throughout the 
magnetosphere, being delayed the interior only the short propagation 
time hydromagnetic effects, the related enhancement precipitation 
would begin with the onset the storm. 

the other hand, particles captured from the strengthened solar wind 
would reach the interior only after the relatively long convection time 
few hours, and would only then provide enhanced rate injection into the 
particle belts. This rate injection would reach maximum when the avail- 
able solar particles and the turbulence reach maximum, presumably about 
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the time the maximum the main phase. would then expect the 
integrated contribution the detectable flux reach its maximum some 
later time still, when the enhanced injection rate returns once again its 
equilibrium value. The anticipated sequence events illustrated Fig.12. 


FLUX (CURVES AND 
INJECTION RATE (CURVE 


TIME (DAYS) 


Fic. The temporal variation the counting rate the outer Van Allen belt (1). 
The curve (2) represents the counting rate particles which formed the belt the time 
the sudden commencement, according the description events given the text, and (3) 
represents the injection rate required produce (1). The variation (3) essentially that 
the storm-time ring current. 


IV. CONCLUDING DISCUSSION 


IV.1. Comparison with Related Theories 

The model have presented appears the first which magneto- 
spheric convection plays important role. Indeed, except for Gold’s emphasis 
and its potential, this convection has been totally overlooked the past. 
nevertheless implicit any model that drives the current system from 
above. Martyn’s (1951) theory prominent representative such models, 
and brief examination contrast our own can instructive. 

Martyn employs the charge separation Chapman and Ferraro his 
driving mechanism, and have already indicated that could serve equally 
for ours. The great extension the earth’s exosphere was not appreciated 
the time Martyn wrote, the convection ionization within played 
part his development. Convection the and regions did, however, and 
Martyn deduced pattern circulation similar our Fig. But his 
derivation the sense flow was opposite ours, and this forced him invoke 
greater mobility for positive than for negative ions order account for 
the sense the current system. Such assumption unjustified theo- 
retically and opposed the correlation bay disturbances with auroral 
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motions and the over-all trend pre-midnight and post-midnight motions 
actually observed. 

The discrepancy between Martyn’s theory and ours can traced the 
fact that Martyn assumed the charge-separation process occur geo- 
centric distance only earth radii, and drive the circulation auroral 
latitudes directly. our model, the viscous-like driving agency acts 
distances greater than earth radii, and operates directly only the 
regions linked zone Fig. 10. The motion across the rest the polar 
cap and auroral latitudes then constitutes return flow, opposite direction 
the flow determined Martyn. have seen, this new pattern fully 
compatible with the observations. 

The theory developed recently Piddington closely related 
ours, that adopts viscous-like interaction produce the geomagnetic 
tail (and indeed, was the first so) and drive the current system. 
both cases, however, Piddington appears employ the viscous interaction 
only deform the field and exospheric ionization from one static condition 
another. therefore fails exploit, even recognize, the important 
consequences that convection can produce. this regard, should noted 
that any distinction between distorted static geomagnetic field and distorted 
convecting geomagnetic field purely semantical distinction, although the 
second these descriptions conjures picture events that more likely 
useful properly applied. the case the ionization, course, the 
distinction between static and convective system perfectly real, and 
are forced believe convection. 

this connection, should comment Piddington’s view the main 
phase geomagnetic storm. understand him correctly, Piddington 
pictures this consequence expansion the low-latitude torus 
the geomagnetic field time diminishing solar-wind pressure, with the 
expansion acting fill gap left field lines that were swept into the geo- 
magnetic tail. would seem that the enhanced viscous interactions must 
maintained the same time, order keep the field lines the tail, 
that some stability the tail must assumed. While agree that this 
process might play some part, would note that the necessary assumptions 
cannot strongly defended and that their failure would certainly limit the 
importance the mechanism. But, what more important, note that 
convection through the interior can return the field lines the tail their 
initial configuration the low-latitude torus, without running counter 
the viscous drag the solar wind (and indeed being aided it). Accordingly, 
have serious doubts that this process could compete effectively against 
the ring current discussed our preceding section. the latter free from 
the difficulties that have been associated with ring currents the past, 
consider likely the most dominant mechanism. 

need not pursue this comparison theories further, except note 
again that mechanisms capable driving the current system from high 
the magnetosphere should lead convective motions the type have 
discussed. Theories auroral and geomagnetic storms that not include 
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such mechanisms may also imply convection, but the circulatory patterns 
would likely different form orientation; their comparison with the 
present model lies beyond our immediate scope. 


IV.2. Summary 

have discussed here mechanism for the generation convective 
motion the outer magnetosphere and have examined very broadly some 
its implications. The latter are good agreement with, and provide forth- 
right explanation of, wide variety features high-latitude geophysical 
phenomena and storm-time variations. are led believe that the con- 
vection system major importance these phenomena, and expect 
provide new basis which theories detail may time based. 
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STUDY THE EFFECT GAMMA RADIATION 
THE MOLECULE 


ABSTRACT 


Modern instrumental techniques were used investigate the effect gamma 
radiation the molecule pentobarbital. The results were correlated with each 
other and with those from biological test. The main effect radiation appeared 
centered about the two substituents the molecule the number 


INTRODUCTION 


preliminary experiment had shown that gas was produced when pento- 
barbital was irradiated the dry state Co™ gammas. Identification the 
many components this gas corroborated the radiation damage mechanism 
suggested electron paramagnetic resonance study pentobarbital 
prepared and irradiated boric acid glass (Cloutier 1961). the purpose 
this report present the experimental findings and correlate them with 
other data from chemical and biological tests gamma-irradiated pento- 
barbital sodium. 

the Samples 

The samples were irradiated 500-curie Transfer Case Irradiator 
obtained from the Commercial Products Division the Atomic Energy 
Canada Limited. The average dose rate was 220,000 Absolute 
dose-yield determinations were not attempted this experiment. Samples 
were simply irradiated level which was known experimentally yield 
concentrations which allowed study the given effect. Exposure times 
varied from few hours long 12} months. 

For the gas analysis experiments, pentobarbital was irradiated sealed 
pyrex container which had been either evacuated filled with nitrogen 
atmospheric pressure. The pyrex containers were designed permit sampling 
the gaseous phase with minimum losses contamination. 


Material 
The free acid and sodium salt pentobarbital were commercial products 


high purity. Levi and Hubley (1956) used samples from the same lot number 
characterize pentobarbital and study its copper-pyridine complex. The 
4,6-di-p-nitrobenzyl (Chatten and Levi 1957; Cloutier and Manson 1959), the 
1,3-dimethyl (Manson and Cloutier 1961), and the 4,6-dixanthyl (Flann and 
Cloutier 1961) derivatives pentobarbital were prepared and characterized 
members this laboratory. Specimens from the same preparations were 
available for this study. 

received June 26, 1961. 
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nembutal 


pentobarbital 
pentobarbitone (B.P.) 
barbituric acid 


—CH 
sodium pentobarbital 
N==C 
3 
pentobarbital pentobarbital 


Fic. Structural formulas pentobarbital and some its derivatives. 


The structural formulae the above compounds are given Fig. 
will noted that the structural formula for the dixanthyl derivative differs 
from that suggested McCutcheon and Plein (1949). recent infrared study 
the reaction barbiturates with xanthydrol (Flann and Cloutier 1961) 
indicated that the two xanthyl groups were attached positions 4,6 the 
malonylurea ring and not positions 1,3 previously assumed. 


INFRARED DATA 
The gas produced during irradiation pentobarbital, exposed months 
gammas, was analyzed with Model infrared spectro- 
photometer equipped with sodium chloride optics. The gas sample was intro- 
duced into evacuated 7.5-cm infrared cell having sodium chloride windows 
and the gas pressure inside the cell was adjusted atmospheric pressure 
with dry nitrogen. The spectrogram Fig. was obtained minutes 
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resolution setting 927. The results reported Table were based cor- 
relations with Sadtler standard gas spectra, the spectra Pierson al. (1956), 
and the spectra from gases obtained from the Matheson Company. 

Two other experiments were made with pentobarbital sodium irradiated 
for weeks and months. The same gaseous components were identified; 
however, the relative proportions were not exactly the same. 


GAS CHROMATOGRAPHIC DATA 


Measured volumes the same radiation-produced gas were also analyzed 
with Model 154-C gas chromatograph. Two types 
Elmer partition columns were used: column diisodecyl phthalate and 
column silica gel. The operating conditions were generally follows: 
helium carrier inlet pressure; 25° 2-meter column; detector voltage, 
volts. 

The results pertaining the irradiation pentobarbital are 
presented Table For identification purposes, the observed elution times 
were compared with those gases generated and obtained from 
either the Matheson Company (CO and the Phillips Hydro- 
carbons Company (ethane, propane, and iso-butane, and iso-pentane). 
The estimates relative concentrations were made comparing the peak 
areas the sample chromatograms with those standard ones obtained with 
known volumes gases. 

The same gaseous components were found following the 3-week irradiation 
pentobarbital sodium. 


ELECTRON PARAMAGNETIC RESONANCE DATA 


The effect radiation pentobarbital and five its derivatives was 
studied room temperature with the method electron paramagnetic 
resonance. The spectrometer used has already been described (Cloutier 1961). 
Irradiation pentobarbital air nitrogen atmosphere, Figs. and 
3D, did not give rise any detectable signal. Pentobarbital sodium, when 
irradiated air nitrogen atmosphere, Figs. and 3H, exhibited 
weak singlet, oersteds wide, which disappeared within minutes after 
irradiation. The dose level had marked effect the intensity, the shape, 
and the rate disappearance the observed signals. 

The microwave absorption curves Figs. 3A, 3E, 3J, and were 
obtained from pentobarbital, pentobarbital sodium, 1-3-dimethyl pento- 
barbital, 4,6-dixanthyl pentobarbital, and 4-6-di-p-nitrobenzyl pentobarbital 
respectively, when prepared and irradiated boric acid glass. This rigid 
glassy dispersive medium has been shown stabilize room temperature 
the free radicals produced radiation many barbituric acids (Cloutier 
1961). The above signals remained unchanged over different experiments 
and for more than years these experiments. seen that the same 
signal, odd-line spectrum, 110-oersteds total width, was obtained from the 
irradiated free acid, Fig. 3A; its sodium salt, Fig. 3E; and its dimethyl deri- 
vative, Fig. interesting note that the first two compounds were 
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Fic. resonance absorption signals from gamma irradiated pentobarbital 
and some its derivatives room temperature. First derivative curves. Markers are 
oersteds apart. Vertical arrows indicate the position D.P.P.H. resonance. Photographs 
and pentobarbital irradiated boric acid glass; photographs and pentobarbital 
irradiated air and nitrogen respectively; photographs and pentobarbital sodium 
irradiated boric acid glass; photographs and pentobarbital sodium irradiated 
air and nitrogen respectively; photograph 1,3-dimethyl pentobarbital irradiated boric 
acid glass; photograph pentobarbital irradiated boric acid glass; photo- 
graph pentobarbital irradiated boric acid glass; photograph 
the signal photograph integrated. 
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crystalline materials while the last one was oil room temperature (Manson 
and Cloutier 1961). The similarity the results with solutes two different 
physical states further experimental proof the versatility and efficiency 
the boric acid glass method preparation for the stabilization room 
temperature the free radicals produced radiation barbituric acids. 
The xanthyl and p-nitrobenzyl derivatives, Figs. and 3K, however, gave 
unresolved triplet about oersteds wide. Qualitatively, the relative 
importance the xanthyl p-nitrobenzyl groups molecular weight basis 
changes the configuration the molecule and would expected modify 
the radiation damage mechanism. These two derivatives were included mainly 
supplement the electron paramagnetic resonance study irradiated pento- 
barbital but the results were too complex for easy interpretation. 

The patterns Figs. and were repeated much higher amplifi- 
cation and are shown Figs. and respectively. The main structure 
the signals was off scale; however, the existence and positions the weaker 
outer components were thus ascertained. The signal Fig. was integrated 
and shown Fig. 3L. The sloping base line was due the slight asymmetry 
the original first derivative curve. 


BIOLOGICAL TEST DATA 


Fifty Wistar rats, weighing 114 grams, were sorted randomly according 
weight into five equal groups. Pentobarbital sodium, irradiated 39, 
98, 150, and 300 Mrad, was injected subcutaneously, one irradiation level per 
group, dose 22.5 mg/kg body weight. The drug was administered 
0.75% solution water. The dose 22.5 mg/kg body weight was selected 
after elaborate preliminary study the dose-response characteristics 
irradiated and non-irradiated pentobarbital sodium with rats. The criteria 
used for determining the onset and the end sleep were based the reported 
methods observation Moir (1937), Carmichael (1938), and Brook and 
Cartland (1944). 

The average sleeping times observed within groups are listed Table 


TABLE 


Average sleeping time rats given dose 
22.5 mg/kg body weight irradiated pento- 
barbital sodium 


Average 


Radiation sleeping 
dose time Log, average 
(Mrad) (min) sleeping time 
3.75185 
3.50856 
150 2.84491 
300 5.34+1.9 1.66771 


along with their calculated standard errors. The log, the average sleeping 
time versus the radiation dose pentobarbital sodium plotted Fig. 
was found that, for dose mg/kg body weight, with rats, linear 
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LOG SLEEPING TIME 


100 200 300 
RADIATION 


Fic. logarithm the average sleeping time rats given dose 22.5 mg/kg 
body weight irradiated pentobarbital sodium. 


regression the log, the average sleeping time (7) radiation dose 
pentobarbital sodium (D) was significant less than 1%, the regression 
equation being: 3.00928—0.00697 D(Mrad) 117.4. 


OTHER TESTS DATA 

The X-ray powder diffraction patterns pentobarbital sodium irradiated 
progressive levels from 300 Mrad were obtained. The spectra showed 
significant changes the positions and relative intensities the lines. 

Pentobarbital sodium, irradiated and 150 Mrad, was assayed non- 
aqueous titration observe any change potency due irradiation. The 
samples were dissolved glacial acetic acid, two drops crystal violet 
were added, and the solutions were titrated with perchloric 
acid dioxane, 0.1 Values 97.8% and 97.2% were obtained from 
unirradiated and radiated pentobarbital sodium respectively. 

infrared differential absorption spectrum pentobarbital sodium 
irradiated 150 Mrad versus the unirradiated salt was obtained. Minor 
differences only were indicated the differential spectrum. These differences 
could not interpreted terms molecular changes irradiated pento- 
barbital sodium. 

water solution 2.68 pentobarbital sodium irradiated 150 Mrad 
was prepared. insoluble fraction was isolated filtration. represented 
3.62% weight the original solute. This water insoluble fraction had 
broad melting point range whose lower limit was 35.0° Its sodium content, 
determined flame photometry, was not significantly different from that 
the water soluble fraction. 
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DISCUSSION 


The electron paramagnetic resonance spectra recorded from irradiated pento- 
barbital, pentobarbital sodium, and 1,3-dimethyl pentobarbital, Figs. 3A, 
3E, and respectively, are odd-line spectra. Their contour characteristics 
indicate hydrogen hyperfine interaction with the electron spin. Further, the 
observed line intensity ratios the above spectra strongly suggest that the 
even number interacting hydrogen atoms are the equivalent the near- 
equivalent type. The spectra Figs. and 3F, obtained much higher 
amplification, show that signal multiplicity high nine must con- 
sidered even though, doing so, the line intensity distribution the signals 
departs seriously from true binomial one. However, believed that the 
observed absorption signals are complex. They may result from the super- 
position signals varying odd multiplicities, centered about the same 
g-value and having the same line spacing (12 oersteds). The relative amounts 
each contributing free radical species would such produce the 
observed line intensity distribution the signals. Restricting the discussion 
the case where single hit per molecule occurs, the various possibilities 
with pentobarbital, taking into account the above experimental data and their 
theoretical implications, are listed Table III. The solid phase free radicals, 


TABLE III 
List the free radicals produced gamma irradiation pentobarbital 


Solid phase Gaseous phase 
free radicals free radicals Type 
CH; 


*H radicals may originate from the destruction of the malonylurea ring structure. 
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trapped the boric acid glass matrix, are all expected give resonance 
absorption signals which have line spacing oersteds (Cloutier 1961). 
Their g-values are also very close that the free electron. Therefore, the 
additive property any given mixture such contributing signals yield 
the observed absorption patterns must considered. The gaseous phase 
free radicals can diffuse and recombine number ways, illustrated 
Table IV. interesting note that Table explains the formation 


TABLE 


Recombination the gaseous phase free radi- 
cals produced gamma irradiation pento- 


barbital 
Recombination 
reaction 
(types are identified 
Table Gas produced 

CH, 


all the gaseous components identified gas chromatography and infrared 
spectrophotometry, Table with the exception and CO. and 
are believed originate from cases where the malonylurea ring structure 
the barbiturate opened the incident radiation. The fragments would 
degenerate rapidly release some any solid residue would 
accounted for the 3.62% water insoluble fraction, already reported for 
irradiated pentobarbital sodium. n-Hexane, listed Table IV, was not 
detected this experiment. However, quite possible that, produced, its 
concentration was not high enough allow positive detection. 

The partial destruction radiation the side chains position 
pentobarbital sodium, outlined Table would account for the decrease 
average sleeping time values function dose level observed the 
biological test with rats. Qualitatively, reduction the complexity the 
substituents position decreases the time during which the barbiturate 
metabolized the body, i.e. the length action the drug. interesting 
note that logarithmic relation between average sleeping time and radiation 
dose pentobarbital sodium was found, would expected with the present 
radiation damage model: reduction average sleeping time directly related 
radiation damage which logarithmically related radiation dose. 

The partial destruction the side chains position would also explain 
the results from the non-aqueous titration irradiated pentobarbital sodium. 
The molecules with modified substituents position would titrated 
barbiturate just the unmodified one: therefore, change potency 
would detected. The fragments corresponding the cases where the 
malonylurea ring opened (the water insoluble fraction) could either 
have been titratable barbiturate possibly have escaped the sensitivity 


the method. 
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CONCLUSION 


chemical bond rupture produced anywhere this molecule, irradiation 
pentobarbital sodium gamma radiation can either disrupt the ring 
structure the side structure position the ring. When comparing the 
results the biological test and the water solubility test pentobarbital 
sodium irradiated 150 Mrad (60% reduction average sleeping time with 
rats opposed only 3.62% water insoluble residue), would seem 
evident that radiation damage the side chain position for one reason 
another preferred process. The relative concentration reported 
Table appears contradict this statement. However, the many 
recombination processes for the gaseous phase free radicals listed Table IV, 
very important possibility has not been mentioned previously, that is, the 
capture slowly diffusing gaseous phase free radicals neighboring solid 
phase free radicals. Such exchange reaction would result change 
physiological activity the irradiated drug but not increase the 
amount alkane gases produced. Therefore, the relative concentrations 
and all the alkane gases, reported Table may not truly 
indicative the relative frequency radiation damage the two molecular 
areas (ring side structure). 

The site bond rupture along the side structure the molecule pento- 
barbital restricted positions where even number nearly equivalent 
hydrogen atoms left; therefore, the complete removal radiation either 
one substituent not expected. immediate consequence that, 
since the physiological activity this type compound known 
associated with the 5,5-disubstitution alkyl, aryl, alicyclic groups (Levi 
1957), radiation, neglecting the less frequent cases where the malonylurea 
ring opened, not expected destroy the physiological activity pento- 
barbital but simply modify it. There definite indication preferred 
positions along the two side chains where bond ruptures may effected. The 
production and ethane has most likely been enhanced play 
diffusion rates and average recombination times the gaseous phase free 


radicals. 
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CHANGES THE DIURNAL HOUR MAXIMUM 
THE COSMIC-RAY INTENSITY! 


ABSTRACT 


The diurnal hour maximum the meson component changed progressively 
Ottawa, Canada, from min min during the period January 
1955 December 1960 while the nucleon component changed from min 
min for the same period. This evidence favors the 22-year cycle 
the diurnal hour maximum that was first suggested Thambyahpillai and 
Elliot, for stations within geomagnetic latitude belt between 58.1° and 
The diurnal hour maximum Churchill changed from min 
min during the period April 1957 December 1960 for the meson 
component and from min min for the nucleon component. 
Although the change was for later hour the indication 22-year cycle 
Churchill not impressive. Resolute the diurnal hour maximum domi- 
nated the varying magnetic masses interplanetary space. shown that 
the anisotropy varies both magnitude and direction depending the con- 
ditions that exist the solar system. 


INTRODUCTION 


The intensity the primary cosmic-ray particles arriving the earth 
not completely isotropic. small anisotropy about 0.3% solar time does 
exist stations not within the Arctic circle and order determine the 
source cause this anisotropy, the solar daily variation has been studied 
extensively many authors.* The possible existence 22-year cycle 
the diurnal hour maximum for the period 1932-1952 similar the 22-year 
cycle the polarity sunspots was demonstrated Thambyahpillai and 
Elliot (1953) for eight stations geomagnetic latitude belt between 58.1° 
and attempt was made Brunberg and Dattner (1954) deter- 
mine the general direction the anisotropy for these stations after correcting 
for the deflection the primary particles the earth’s dipole field. They 
concluded that the daily variation caused excess particles arriving 
from the 18-hour direction and due the tangential component the solar 
wind. attempt was made explain the hour maximum obtained 
Manchester, which does not agree with this conclusion. the cycle the 
diurnal hour maximum repeated the present sunspot polarity cycle 
station other than Manchester then the above explanation Brunberg 
and Dattner the diurnal hour maximum not adequate. Ottawa and 
Manchester are the same geomagnetic latitude and hence comparable. 

The present paper deals with the results obtained Ottawa, Churchill, 
and Resolute. Ottawa (56.8° N.) within the latitude belt mentioned above, 
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station. Resolute the cutoff magnetic rigidity effectively zero and the 
magnetic lines force move out into space beyond earth radii and become 
entangled with whatever magnetic fields exist interplanetary space. Thus 
anisotropy always exists the low energy portion the primary particle 
spectrum will apparent Resolute time low magnetic turbulence 
interplanetary space. 

Clark al. (1961) have shown that primary particles energies 
are isotropic. anisotropy that due the sun will due particles 
much lower energies and change the diurnal hour maximum can 
expected with changes the sun’s activity. Also, the solar dipolar field which 
was considered fixed polarity now known reverse its polarity (Babcock 
1961). important, therefore, determine whether the 22-year cycle 
the diurnal hour maximum genuine and, so, this will add weight 
the hypothesis that the daily variation extraterrestrial. 

will shown that there strong evidence Ottawa that the beginning 
the 22-year cycle the diurnal hour maximum repeating itself the 
present sunspot polarity cycle but that the effect Churchill not impressive. 
There indication 22-year cycle Resolute. There evidence, 
however, that the diurnal hour maximum Resolute changes with changes 
the planetary magnetic index 


EXPERIMENTAL RESULTS 


The cosmic-ray intensity measured neutron monitors influenced 
the changes the atmospheric pressure, whereas the meson component 
influenced the distribution the mass air above the recording instru- 
ment. This mass distribution will depend upon the time year and hence 
12-month moving averages are used this effect will negligible. 

Twelve-month moving averages the pressure corrected data from cubical 
telescopes and neutron monitors Ottawa, Churchill, and Resolute were 
analyzed harmonically. The amplitude and phase the first harmonics 
the meson component are shown Fig. can clearly seen that the 
diurnal hour maximum Ottawa becomes progressively later the solar 
cycle progresses but does not follow the activity the sun. The hour 
maximum Churchill changes minutes while for the same period the 
change hours minutes Ottawa. The diurnal hour maximum 
Resolute appears move through cycle. 

The amplitude and phase the nucleon component for the three stations 
are shown Fig. 1955, the neutron monitor Ottawa contained two 
counters. The design the pile was changed and the new pile contained three 
larger counters (Fenton, Fenton, and Rose 1958). The diurnal hour maximum 
for 1955 given but not included the 12-month moving averages which 
were started with the data from the new pile. can accepted that the 
diurnal hour maximum the nucleon component followed that the 
meson component between 1955 and 1956 since the hour maximum changed 
hour minutes for the nucleon component and hour minutes 
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Fic. moving averages the amplitude and hour maximum the 
first harmonic the meson component and the activity the sun measured the solar 
10.7-cm radio flux, supplied Covington. 


for the meson component and both cases the change was towards later 
hour. The same apparent cycle present the nucleon component Resolute. 


DISCUSSION 


The daily variation the cosmic-ray intensity has been considered 
due anisotropy the cosmic rays and the rotation the earth. Since 
the earth possesses magnetic field and this field only considered then the 
diurnal hour maximum will vary one moves from the equator the 
pole, due the deflection the charged primaries the earth’s field. Hence, 
the anisotropy, the magnetic field, and rotation the earth are the only 
factors influencing the diurnal hour maximum then the hour maximum 
should remain constant any particular station provided there change 
the magnetic rigidity spectrum the particles and significant change 
the earth’s magnetic field. That the cosmic-ray intensity changes with solar 
activity has been amply demonstrated Forbush (1954) for two solar 
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Fic. Twelve-month moving averages the amplitude and hour maximum the 
first harmonic the nucleon component and the activity the sun measured the 
solar 10.7-cm radio flux, supplied Covington. 


cycles. The earth’s field, though disturbed, can considered constant for the 
short period under study. The intensity the low magnetic rigidity particles 
reduced with increase solar activity and, therefore, the diurnal hour 
maximum should either remain constant change earlier hours. 
Thambyahpillai and Elliot (1953) have shown that the change the 
diurnal hour maximum for the meson component has approximate 
22-year cycle similar the 22-year cycle the polarity sunspots. For 
approximately eight years, 1939 1946, the hour maximum was 
15+0.5 and either end the cycle about min. The 
hour maximum for the meson component Ottawa changed from 
min 14hr min from January 1955 December 1960. The change 
the diurnal hour maximum was progressive with significant abrupt 
change. The nucleon component changed from min min. 
The difference the hour maximum between the two components 
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negligible for 1960 but very significant for 1955. Both components reached 
the same diurnal hour maximum approximately hours time 
decreasing solar activity though the intensity the nucleon component has 
not increased significantly above the low the time maximum solar 
activity. Since considered that the intensity low magnetic rigidity 
particles decreased thus causing the decrease the general intensity, the 
diurnal hour maximum shouid not progress later hours without corre- 
sponding increase intensity. Hence obvious that the anisotropy not 
directly influenced the earth’s magnetic field Ottawa. The direction 
the primary particles must influenced fields force external the 
earth’s magnetic field and any anisotropy would due the resultant 
all the fields. This resultant field will different for different parts the 
world. 

The diurnal hour maximum Resolute the nucleon and meson com- 
ponents differ continuously approximately two hours. The apparent cycle 
which, for this epoch, appears have period almost two years present 
both components and the planetary magnetic index with the minima 
and maxima preceding three months those the cosmic-ray diurnal hour 
maximum, Fig. Since the cycle the hour maximum does not appear 
either Ottawa Churchill must concluded that the diurnal hour 
maximum Resolute more sensitive the interplanetary magnetic dis- 
turbances measured the planetary magnetic index 

differentiate between the mechanisms that influence the daily variation 
important determine (1) whether not the diurnal hour maximum 
remains the midafternoon stations within geomagnetic latitude belt 
between approximately 60° and 60°S., the next minimum solar 
activity, (2) whether not the anisotropy stations near the geomagnetic 
poles controlled continuously the magnetic state interplanetary space. 
these conditions prevail then the anisotropy caused prime mechanism 
with two different effects (i) anisotropy related the 22-year cycle 
the general polarity the sun and observed stations within geomagnetic 
latitude belt between 60° and 60°S., (ii) anisotropy polar stations 
controlled some magnetohydrodynamic effect interplanetary space. 

determine the diurnal hour maximum Resolute time low 
interplanetary disturbances, the data from second telescope stationed there 
were analyzed. This telescope similar the one used aboard the Canadian 
icebreaker (Rose and Katzman 1955) measuring cosmic rays 
during cruise through the Northwest Passage. The amplitude and phase 
the first harmonic shown Fig. together with the data from the I.G.Y. 
cubical telescope. The gap the data from the smaller 
due faulty operation the equipment and the fact that the 12-month 
moving averages were started from March 1956. This telescope was shut 
down early January 1960 when was found that trouble had developed 
sometime earlier. The diurnal hour maximum for the year November 1954 
October 1955 that was observed, Resolute, when was very low, was 
02hr with amplitude only place very little 
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Fic. Twelve-month moving averages the diurnal hour maximum the meson 
and nucleon components Resolute and the monthly mean the 12-month moving averages 
the planetary magnetic index shifted months. 


reliance phase angles when the amplitude low. The amplitudes obtained 
from the data both telescopes are the same for the period when both were 
operating; permissible, therefore, accept the result the smaller 
telescope for the earlier period. Also the hours maximum follow each other 
fairly well except for the last months when sporadic trouble developed the 
smaller telescope. Hence considered that, time low interplanetary 
magnetic disturbance, anisotropy does exist Resolute the meson 
component would the order 0.02+.01% hr. With increasing 
interplanetary magnetic disturbance measured the intensity the 
anisotropy increases Resolute and the hour maximum changes. 

The amplitude the diurnal variation the meson component Ottawa 
increased 100% from 1955 1958 and then decreased approximately 
the 1955 value (Table while the 1960 value the nucleon component 
50% above the 1955 value. Churchill the changes the amplitudes are 
not very significant. Resolute, though the amplitude the first harmonic 
the nucleon component remains constant, the amplitude the meson 
component increased 100% between 1957 and 1960. the magnetic 
rigidity spectrum simply depressed during the time solar maximum 
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Fic. moving averages the amplitude and diurnal hour maximum 
the meson component obtained from the data the and 
telescopes Resolute. 


activity, discussed Simpson (1960), then the amplitude the diurnal 
variation should remain constant. This approximately true for the nucleon 
component but not for the meson component. Since particles all energies 
can arrive Resolute and the air mass the only large absorber, then the 
amplitude the nucleon component should higher than the meson com- 
ponent, due the decay the meson. This for 1957 and 1958 but not 
for 1959 and 1960 when the amplitudes are equal. 

accepted that Resolute the diurnal hour maximum modulated 
the interplanetary magnetic masses, permissible obtain mean value 
for the period October 1956 December 1960. The mean diurnal hour 
maximum for the meson component +20 min and min 
+24 min for the nucleon component. Churchill, the mean hour maximum 
for the meson component min and min for the nucleon 
component. Ottawa, the hour maximum for both components becomes 
progressively later and for 1960 the values are min and min 
for the meson and nucleon components respectively. obvious from the 
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above discussion that the direction the anisotropy cannot obtained 
considering only the deflection the primary particles caused the earth’s 
field. The anisotropy varies both magnitude and direction depending 
the conditions the solar system. more information about these conditions 
becomes known, clearer picture will emerge with regard the cause 
the observed anisotropy. 


CONCLUSIONS 


The diurnal variation the cosmic-ray intensity Ottawa, Churchill, 
and Resolute not consistent with anisotropy that fixed space 
time. The hour maximum the diurnal variation Ottawa becomes 
progressively later, which indication that the diurnal variation follows 
the 22-year cycle the polarity sunspots was first demonstrated 
Thambyahpillai and Elliot. The anisotropy the meson component Reso- 
lute questionable time low planetary magnetic index but 
definite for high values when the direction the anisotropy becomes 
dominated the interplanetary magnetic turbulence. 
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HIGH-FREQUENCY DIFFRACTION SPHERE! 


ABSTRACT 


systematic procedure given for the determination the asymptotic 
series directly from the Helmholtz equation and the boundary conditions for the 
field the shadow region sphere illuminated plane wave high fre- 
quencies. The first two terms the series for the shadow region, including the 
regions near the axial caustic and the boundary layer near the surface the 
sphere, are explicitly evaluated. The present procedure can generalized 
any smooth convex three-dimensional object with rotational symmetry, illumi- 
nated plane wave the direction the axis rotational symmetry. 


INTRODUCTION 


recent times the problem high-frequency diffraction obstacles 
shape’’ whose geometry does not fit into simple co-ordinate 
system has received considerable attention. Keller (Levy and Keller 1957) 
has treated the problem diffraction convex smooth object using geo- 
metrical theory diffraction that generalization the geometrical optics 
the reflected and refracted rays through the introduction additional rays 
called the the rays. The method Keller provides 
one more term various asymptotic developments addition that given 
geometrical optics. But, without further elaboration, information about 
the higher-order terms can obtained the above method. recent 
paper (Wu 1957), new method has been proposed which 
possible obtain higher-order terms the asymptotic series for the field 
the shadow region circular cylinder through the use the boundary 
layer theory. This procedure was later used and Seshadri (1958) 
solve the problem diffraction plane waves infinite cylinder whose 
cross section smooth convex curve. Together with the method Keller 
al. (1956), which works satisfactorily for the illuminated region, the problem 
diffraction infinite convex cylinder can considered well 
understood. The next problem line for investigation the diffraction 
convex three-dimensional obstacle with rotational symmetry, illuminated 
plane wave the direction the axis rotational symmetry. Again, 
the illuminated region the method Keller can used satisfactorily. 
the shadow region, there appears additional feature which not present 
the case cylinder; that is, axial caustic appears the shadow region 
along the axis rotational symmetry the body. this report the simplest 
case this class obstacles, the sphere, solved the sense that systematic 
procedure given for the determination the asymptotic series for the field 
the shadow region, including the regions near the axial caustic and the 
surface the sphere. The generalization the case arbitrary convex 
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body with rotational symmetry will given subsequent report. Although 
the explicit evaluation for the latter case considerably more complicated, 
there are new basic difficulties principle. 


STATEMENT THE PROBLEM AND THE INTERIOR 
SOLUTION THE SHADOW REGION 
The problem studied follows. For the incident field 
(2.1) exp(—ikz), 


desired determine the asymptotic behavior, tends infinity, 
the total field which satisfies the time-reduced Helmholtz equation 


and the Dirichlet boundary condition 
(2.3) 


sufficient consider the behavior the positive shadow region, 
the sense geometrical optics, for given value 

The variables used the study the field the shadow region, 
except for the region near the axial caustic, are and defined for given 
follows (Fig. 1): 


(u,t) 


CONSTANT 


(0,0) 


INCIDENT 
PLANE WAVE 


Fic. The (u,t) co-ordinate system for the shadow region. 


~ 
CONSTANT 
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2.4 

The origin the co-ordinate system chosen the point intersection 
the positive shadow boundary with the boundary curve the sphere for 
the given 

the present co-ordinates, the wave equation takes the following form: 


t—cot 


The general solution split into slowly varying amplitude factor 
and rapidly oscillating exponential factor where the form 


determined from the partial wave expansion. 
From the commutation relations 


and 


follows that 


From (2.6) and (2.8), the following form supposed for the wave function 
the shadow region: 


(2.9) 
n=0 
When (2.8) applied (2.9), the following results are obtained: 


(2.11) 


where the constants integration and depend 


@ 
q 
i 
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Equations (2.10) and (2.11) and the F,,’s determined continuing the 
iterative procedure are not valid the vicinities or/and cot 
corresponding the regions near the surface the sphere and the axial 
caustic. Consider first the region near the axial caustic. 


THE FIELD NEAR THE AXIAL CAUSTIC 


view the fact that the co-ordinate system, the neighbor- 
hood the axial caustic, taken along the direction which the field varies 
most rapidly, the following change co-ordinate suggested (Fig. 2). 


the co-ordinates, the neighborhood the axial caustic, all the 
rapid variations and its partial derivatives are now only. 


PLANE WAVE 


The co-ordinate system for the axial caustic region. 


€=0 
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With (3.1), (2.8) becomes 


accordance with the spirit the boundary layer theory (Carrier 1953) 
examination (3.2) indicates that the change variable 


(3.3) 

gives the proper co-ordinate stretching. With (3.3), (3.2) becomes 
where 


1+cot? 


se = 


view (2.9), the following form supposed for the wave function the 
region near the axial caustic 


n=0 
Therefore the first three partial differential equations the caustic region are 


(3.6) 
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and the boundary conditions are that all the must bounded for 
and that must blend with the interior solution, the regions where 
both solutions predict the same behavior From (3.9) and the boundary 
conditions, follows that 


sin 
cos 
From (3.10), (3.12), and the appropriate boundary conditions, follows that 


sin 


Jo(n 


(3.13) 


This iterative procedure can continued give all the higher-order terms. 
The remain determined. 


THE FIELD THE BOUNDARY LAYER NEAR THE SURFACE 
THE SPHERE; CONCLUDING REMARKS 


The field the region near the boundary curve the sphere and the various 
constants integration, (2.10) and (2.11) can again determined 
through the use the boundary layer theory and the Dirichlet boundary 
condition. Let 


(4.1) 


Then rearrangement terms the right-hand side (2.8) gives 


2/3 2 20 1 3? } 


The leading term (4.2) can simplified letting 


(4.3) 
Then follows that 
where 
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Let 

where 

n=0 


Once known, and hence and can determined the asymp- 
totic expansion This will illustrated for the case 

The application (4.4) (4.8) and (4.9) gives the first three partial 
differential equations the boundary layer: 


The outgoing-wave solution (4.10) the following Airy function. 

The application the Dirichlet boundary condition (4.13) gives 


From (4.11) and the boundary condition 


follows that 


cot 


Through integration parts, and since the integral the square the Airy 
integral not zero, follows that 


(4.15) Cro(u) 2cot 
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Therefore 


Thus the leading term the asymptotic development for the field the 
boundary layer near the sphere completely determined. 

Before proceeding the determination now interesting 
determine and complete the leading term the interior solution 
and the solution near the axial caustic. From (4.17), using the asymptotic 
representation for the Airy integral and comparison with (2.10), the 
following result obtained. 


From (4.18) and (2.10), the following interior solution obtained. 


From (4.18) and (3.12), the following solution near the axial caustic region 
obtained. 


Two remarks are now order. First, the functional forms and 
given (4.19) and (4.20) are agreement with results either the partial 
wave expansion Keller’s (1957) geometrical theory diffraction. Secondly, 
the fact that the phase decreases ray passes through focus 
placed evidence the factor (cos (4.17), since cos changes 
For the second term the asymptotic series for the field the shadow 
region, only necessary determine From (4.12), (4.17), (4.18), and 
the boundary condition 


follows that 


where 


0 


| 
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The integral summable the sense Abel and hence number (Wu 
and Seshadri 1958). The solution (4.21) gives 


where 

= 2 
(4.25) 
Therefore 


The second term the asymptotic series for the field the shadow region 
thus completely determined. 

Several remarks should made connection with the results obtained. 
Two normalization constants and are still undetermined, although these 
are not prime importance here. This due the fact that the magnitude 
the incident wave has not been used. They can determined through 
comparison with partial wave expansion any one point the shadow 
region. 

Secondly, the results are not valid the neighborhood the foci 
cot 0). The situation here more complicated, since the boundary 
layer theory does not apply. The size the focus region centered 
around any focus. Hence goes zero tends infinity. 
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UNCHARGED CONDUCTING TOROIDAL RING 
UNIFORM ELECTRIC 


ABSTRACT 


Mathematical expressions for the potential function uncharged conduct- 
ing toroidal ring placed uniform electric field are derived and expressed 
terms toroidal functions. Some numerical results were calculated the 
IBM 650 computer the University Toronto and are included the present 
paper. verify the calculated results, systematic study electrolytic 
tank was undertaken. was found that the theoretical calculations agreed well 
with the experimental results. 


INTRODUCTION 


The potential functions uncharged conducting ring uniform electric 
field are usually calculated approximately, assuming that the diameter the 
conductor small compared with that the ring. the present paper the 
exact solutions the problem have been derived means toroidal func- 
tions. Two cases have been discussed, namely, (i) uniform field electric 
force parallel the axis the ring and (ii) uniform electric field perpen- 
dicular the axis the ring. Therefore, the principles superposition 
the general potential functions uncharged conducting ring placed 
uniform electric field any direction may derived. Some numerical results 
calculated IBM 650 computer the University Toronto are presented 
order illustrate the derived solutions. 


UNIFORM ELECTRIC FIELD PARALLEL THE AXIS TORE 
Consider conducting ring radius shown Fig. The ring 
zero potential and situated field that had uniform static field intensity 
before the ring was introduced. That is, the original potential distribution 


given 
(1) 


This must also the distribution for large after the ring introduced. The 
problem, then, find the potential everywhere outside the ring, the 
potential and field being zero inside the tore. 

The potential function electrostatic field system always satisfies 
equation. The general solution Laplace’s equation 
co-ordinates given Carter and Loh (1958) follows: 
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where and are arbitrary constants, 
and are known the first and second 
kinds respectively. 


Fic. Cross section conducting ring. 


The toroidal co-ordinates are related the well-known cylindrical 
polar co-ordinates the following expressions: 


cosh 
(3) 
cosh 


The relations between and are 
arc sech (r/R), 


When the system symmetrical with the central plane, the present 
case, therefore the general potential function becomes 


when approaches infinity. thus get 


| 
| 
n=0 
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and 

sin 


Therefore the solution for the potential function 


The any point the tore 


(8) 


where the scale factor the toroidal co-ordinate system, 


a n=0 


But 
Thus 


UNIFORM ELECTRIC FIELD PERPENDICULAR THE AXIS 
THE TORE 
The general potential function derived (7) longer valid for the 
present case where the electric field perpendicular the axis the ring, 
since this three-dimensional problem. However, the general solution 
Laplace’s equation for the toroidal co-ordinate system (u, has been 
derived the author (1959) follows: 


where dm, Amn, ANd are arbitrary constants. When large, 
the field undisturbed and must equal 


sinh cos 
cosh 


3 
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For the potential vary (12) for large values where the field 
undisturbed the presence the ring, the coefficients and 
must equal zero and also for all values except Thus 
the potential function becomes 


n=0 


where new constant. 
Introducing the boundary condition that 


(14) 
have from (13) that 


1 
cosh 


order evaluate the arbitrary constant multiply eq. (16) 
cos and integrate from 27. Equation (16) yields 


cos 


The integral may evaluated readily using the following identity: 


(cosh 


(18) 


tr 
e 


where are the first derivative Substituting (18) into (17), 
get 

Hence the general expression for the electric potential any point 
outside the ring having uniform electric field perpendicular the 
axis the ring given 


cosh 
44/26, a 
1 
before, the densities charge induced the ring can be, necessary, 
readily evaluated. 


ELECTROLYTIC TANK INVESTIGATION 
order verify the calculated results electrolytic tank investigation 
was carried out. toroidal ring axially symmetric body wedge- 


“Ve 
n 
sinh 
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shaped tank was used with the apex the wedge representing the axis 
the toroidal ring. brass disk electrode was placed such that its distance 
from the wedge apex was equal the toroidal radius The probe used was 
0.010 in. diameter platinum wire. oscilloscope was used for null indica- 
tion. 

preliminary test was carried out determine the accuracy obtainable 
with the experimental setup. The experimental equipotential lines for 
electrode system representing concentric cylinders were compared with the 
known exact solution. was concluded that the experimental error was 
within 1%. 

CONCLUSIONS 

The potential function derived eq. (7) for the electric field parallel 
the axis the ring has been calculated, using the first four terms the 
series, for two cases with different ratio the electronic computer 
IBM 650 the University Toronto. Figures and show the equipotential 


oF ROTATION 


Axis 


RING 


Fic. Uncharged conducting ring originally uniform field (R/r 5). 


AXIS OF ROTATION 


Fic. Uncharged conducting ring originally uniform field (R/r 10). 
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lines near the ring which placed uniform electric field parallel its 
axis. found that the potential field distance nearly undis- 
turbed the presence the ring and independent the ratio R/r. How- 
ever, the field near the ring dependent the ratio R/r. The calculated 
values agreed closely with experiment (1%). 

Though the solutions for the potential function have been obtained separately 
for two cases, namely (i) the field parallel the axis the ring and (ii) the 
field perpendicular the axis the ring, the general potential function for 
the ring placed uniform electric field any direction may derived 
means the principle superposition. 
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SUBSTRUCTURES AND DISLOCATIONS PRODUCED DURING 
SOLIDIFICATION ALUMINUM! 


HERALDO BILONI 


ABSTRACT 


The substructure and the dislocation arrangement produced during poly- 
crystalline solidification aluminum has been studied using the decant- 
ing method for the interruption the solidification. Etching techniques 
reveal individual dislocations and the method thin epitaxial layers aluminum 
oxide reveal changes impurity concentration have been employed simul- 
taneously. The development the structure towards the interior the solid 
was studied successive electrolytic polishings. this way was possible 
show that the ‘‘cell formed dislocations originates areas 
where impurities has occurred. Towards the interior 
the solid dislocation rearrangement into substructure 
which closely related the cellular substructure occurs. the same time 
substructure occurs, arising from the original grain 
boundaries. The results are interpreted the light recent theories solidi- 
fication. 

INTRODUCTION 

Chalmers and his collaborators (Rutter and Chalmers 1953; Tiller and 
Rutter 1956) have described the formation substructures during solidifi- 
cation metals which are closely linked the impurity content (corrugated 
and dendritic substructures) and have explained their origin result 
(Winegard and Chalmers 1954). 

the other hand, Teghtsoonian and Chalmers (1951) have pointed out 
that macromosaic substructure formed monocrystals during solidifi- 
cation which composed crystalline zones which are slightly disoriented 
with respect each other and which have great stability even after prolonged 
annealing. These two characteristics indicate that the boundaries macro- 
mosaics are formed result the alignment dislocations, happens 
the case boundaries subgrains polygonized substructures. 

The theories which have been put forward most frequently explain the 
formation dislocations during the process solidification are: 

(a) Nucleation and ulterior collapse vacancies (Frank 1955). Schoeck 
and Tiller (1960) have proved, however, that this mechanism does not suffice 
itself explain the formation macromosaic striated substructures. 

Microsegregation solute. Tiller (1958), developing ideas originally 
expressed Friedel (1956), supposes that the stages the advancing 
interface (Burton 1957) and the borders the cell sub- 
structure, the microsegregations solute produce variations the crystal- 
line lattice parameter. These distortions are eliminated through adequate 
formation and distribution dislocations. During the process cooling they 
rearrange themselves minimize the energy, thus forming the macromosaic 
substructure. 
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the purpose the present work study the relation existing between 
cell substructure and the formation dislocations, and their subsequent 
rearrangement during the process cooling. The method used com- 
bination the technique etch pits and method using thin epitaxial 
(Biloni and Destaillats 1959), which were deposited anodic oxidation. 


EXPERIMENTAL METHODS 

The samples were melted alumina crucible and cast into nuclear 
graphite mold square section. Prior completion solidification, the 
remaining liquid was drained, exposing the interface (Walton 
and Chalmers 1959). Viewed under the microscope and without previous 
preparation the interface presents appearance similar that the works 
Chalmers and his collaborators, that is, cells whose borders are 
lower level than their centers shown Fig. 

The test tubes were electrolytically polished with Jacquet reagent (1956), 
composed 333 perchloric acid and 666 acetic anhidride, with voltage 
and electrolytic refrigerated vat. 

For the production etch pits, the reagent recommended 
Marchin, and Lacombe (1959) was adopted; the composition used was the 
following: 49% 1.39); 10% 1.185); butylcellosolve, 
40%; AuCl, p.p.m.; HF, 1%. 

The method preparation the reagent and the procedures followed 
were those recommended the above-mentioned authors. With this method 
both micro and geometrical etch pits are developed. This fact, however, does 
not affect the results. The reagent used only once and the period attack 
about one hour. Nevertheless, careful control the time necessary 
each case, since the presence impurities the bath has substantial 
effects upon the time attack (Wyon al. 1959). 

order obtain conclusions from the data furnished the techniques 
corrosion figures, must kept mind that the case aluminum 
(Lacombe and Wyon 1954) (a) etch pits reveal the atmospheres formed 
certain solutes around dislocations; for 99.99% aluminum, the works 
Wyon and Lacombe (1954) prove semiquantitative manner that the 
responsible atom the atom. 

The detection segregation substructures was carried out applying 
the method thin epitaxial layers that show interference colors under oblique 


light (Lacombe and Mouflard 1953). the case aluminum, anodic oxidation 


interface resulting from the interruption solidification, when draining 
the remaining liquid. Cellular substructure with borders enriched solute. 

interface slightly polished with Jacquet reagent and anodically oxidated 
order form thin epitaxial layers. Cellular substructure may detected the different 
color the borders compared with the centers, result the concentration solute. 
may seen that one cell shares two more grains, because the migration the grain 
boundary during the process cooling. X80. 

Fic. Same zone shown Fig. situated microns towards the inner core the 
grain. may seen that the epitaxial layers only indicate the persistence the vertices 
the cells. The continuous line corresponds the so-called macrosegregation, survival 
the original size the grain. X80. 

Interface slightly polished with Jacquet reagent and attacked with fluorated aqua 
regia. Formation etch pits along borders cellular substructure. 


a 


PLATE 


Fics. 


PLATE 


BILONI: SOLIDIFICATION ALUMINUM 1503 


time varies between and seconds, but since additive, possible 
control perfectly well order obtain maximum contrast coloration. 
order study the relation between etch pits and segregation sub- 
structures, the sample attacked with fluorated aqua regia reactive was sub- 
jected anodic oxidation under conditions operation identical with those 


described. 
All micrographs were taken dark field illumination, thus contributing 


the better detection corrosion microfigures and remarkable improve- 
ment the formation interference colors. 


RESULTS THE EXPERIMENT 


Segregation Substructures 
have limited our study substructures cellular character, also 


called Smialowski substructures, the interfaces which present the aspect 


shown Fig. 
Electrolytic polishing few microns, followed anodic oxidation, makes 


possible determine, the difference the interference colors, the con- 
centration solute along the border the cells (Biloni and Destaillats 1959). 
The method epitaxial layers indicates differences orientation the 
crystalline structure, since the thickness the layers determined this 
orientation. This may appreciated Fig. also the superposition 
the cellular network the Smialowski substructure, which may detected 
the different coloration compared with the matrix. 

The explanation for the detection segregated areas means the 
method thin epitaxial layers not yet established, and the possible mech- 
anisms are two: 

(a) Variation the lattice parameters the matrix due impurities 
can change the kinetics anodic oxidation and lead the formation 
different thicknesses the layers (Schoeck, private communication). 


Fic. Zone the interface prepared like the one shown the preceding figure. Formation 
etch pits along the borders the cells. may observed that the response the chemical 
attack varies each grain. the one shown the upper part, the attack revealed segregated 
areas besides producing micro etch pits. X80. 

Fic. situated 200 microns from the interface, attacked with Wyon and Lacombe 
reagent for several minutes. Concentration etch pits the vertices the cells and zones 
corresponding macrosegregation. X80. 

Fic. Zone situated several hundred microns away from the interface, electrolytically 
polished, attacked with fluorated aqua regia reactive, and then anodically oxidated. Grain 
presents strong, heavy concentration etch pits the vertices the cells and alignments 
microfigures that delimit the macromosaic substructure, closely linked with the borders 
the cells detected epitaxial layers. Grain presents concentration micro etch pits 
the vertices the cells and macromosaic substructure delimited alignment micro 
etch pits, which present different sizes though always passing through the vertices the 
cells. The concentration micro pits appearing the center the picture corresponds 
macrosegregated region comprising grains and Grains and have suffered 
practically attack from the fluorated aqua regia reactive and the epitaxial layers reveal 
different stages detection the cell substructure. 

Fic. close the former, but situated 100 microns towards the interior the 
grains. The sample has been electrolytically polished and attacked with fluorated aqua regia 
reagent. Concentration etch pits the vertices the cells, which has produced over- 
attack. Micro etch pits aligned macromosaic substructure, generally passing through cell 
vertices. X80. 


. 
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(b) Different local composition can cause the difference thickness the 
layers, which this case would result from chemical reactivity different 
from that the matrix. 

known that each grain grows with type cell substructure the 
morphology which depends upon crystalline orientation (Tiller and Rutter 
1956). Fig. possible observe the presence two more grains 
that share the same substructure. course, this case, when the border 
cell passes across grain boundary, the color the two sides the 
latter different: the action segregation added that different 
crystalline orientation. These facts are explained take into account the 
possibility migration the grain boundary during the period cooling 
after solidification. 

the electrolytic polishing continued towards the interior the grains, 
the appearance the segregation the borders the cells changes after 
‘several 10-micron polishings. While the cellular substructure persists some 
grains, others the borders the cells disappear and only the vertices 
remain. There wide scale variation between these two cases, and 
there always intensification color the cell vertices. This inter- 
preted indication the fact that when there has been sufficient homo- 
genization, the detection segregations through thin epitaxial layers depends 
large extent upon the orientation the surface each grain. The per- 
sistence the vertices the cells explained there being regions where 
three hexagons meet, i.e., segregation greater than along the sides. 

Figure illustrates the existence substructure consisting continuous 
segregated zones, connected with the grain boundaries the substructure 
and forming net within which there whole number cells. the 
other hand, analysis Fig. makes possible verify that the actual 
boundaries the grain have attained their position through migration during 
the process cooling. 

conclude that the borders the macrosegregation substructure corre- 
spond the limits the cells which originally formed the boundary the 
grain. 

During the process cooling, the mobility the grain boundary has been 
large enough produce migration the present positions, but segregation 
has remained proof their original position. 


Relation between Etch Pits and Segregation Substructures 

The purpose the superposition the techniques etch pits and thin 
epitaxial layers study the relation between segregation substructures 
and dislocations existing the metal the time solidification, taking 


Same zone shown Fig. Brief attack with fluorated aqua regia followed 
anodic oxidation. Geometrical etch pits concentrated along borders macrosegregation 
substructure and some cell vertices. 

The same zone that the preceding picture, but subjected anodic oxidation 
following attack etch pits. Owing the interference colors may seen that while the 
vertices the cells are perfectly clear (because either their color the concentration 
etch pits which have produced overattack), the borders are much more diffuse. Besides, 
may observed that the alignments etch pits generally pass through the vertices the 
cells but not necessarily through their borders. 


III 


| 
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into special account that the etch pits detect dislocation-associated impurities. 
Our experimental results correspond the successive application both 
techniques the utilization etch pits only, and this point clearly 
indicated each case. this work have studied the distribution etch 
pits starting from the interface and their relation with segregated 
areas. 

(a) Interface 

The interface, slightly polished with Jacquet reagent and attacked with 
fluorated aqua regia reagent, presents strong concentration etch pits 
the borders the cells the Smialowski substructure (Figs. and 5). 
may observed the figures that the response the chemical attack 
the result the crystalline orientation the grains, and after prolonged 
attack even the substructure segregation was revealed one them 
(Fig. 5). 

(b) Relation between Segregated Zones, and Etch Pits According 

Distance from the Interface 

Just the aspect segregated zones changes according the distance 
from the interface, the same happens. with the distribution etch pits. 
several tens microns observe strong concentration the macro- 
segregated zones and the vertices the cells (Figs. and 7). more 
prolonged attack produces such concentration points attack the 
vertices the cells that they produce distinct hole shown Figs. 
and 10. 

The distribution etch pits varies from one grain another according 
their crystalline orientation. observe some, starting from distance 
about 100 microns from the interface, marked macromosaic substructure 
revealed the alignment etch pits. This substructure closely linked 
with the cellular substructure (Hulme 1954; Takaka al. 1955), because the 
etch pits usually pass through the vertices the cells and keep certain 
relation with their walls, though not following them closely (Figs. and 
10). The size the macromosaic substructure varies from that corresponding 
one cell that various cells. 

accordance with the concepts set forth Tiller (1958), the foregoing 
results may interpreted follows: 

(1) the zone the interface and the immediate surrounding zone, 
etch pits make possible detect the associations between impurities 
segregated the borders the cell substructure and dislocations formed 
the mechanism microsegregation solute. 

(2) The homogenization produced during the process cooling liberates 
great number dislocations, especially those formed the stages growth 
the interface and along the sides the hexagonal substructure, which 
are eliminated others contrary sign reorder themselves following 
the mechanism climbing suggested Tiller (1958). Our results are 
agreement with the recent work Damiano and Tint (1961). 

strongly segregated zones (vertices cells and zones macrosegrega- 
tion) neither homogenization nor migration the grain boundary liberates 
anchored dislocations, for which reason attack with fluorated aqua regia will 
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reveal almost every case the existence etch pits those zones con- 
centration. 

Liberated dislocations, which have not been eliminated others 
contrary sign, will tend form minimum energy distributions (macromosaic 
substructure) and seems logical suppose that they have connection with 
the existence segregated zones, because the deformation energy edge 
dislocations may diminished solute atoms, different size from 
that the matrix, segregated the dislocations (Cottrel 1948). The same 
mechanism may applied the case alignments dislocations forming 
the substructure the macromosaic (Atwater and Chalmers 1957). 


IV. CONCLUSIONS 

Experiments with foundry polycrystals aluminum, 
solidification has been suspended the sudden draining the remaining 
liquid, furnish ample information about segregation substructures apply 
the method epitaxial lavers producing interference colors under oblique 
light. 

means this technique has been possible determine the variation 
segregation starting from the interface and proceeding towards the inner 
part the grains. Detection, means thin epitaxial layers, the per- 
sistence segregation the vertices the cells and its disappearance the 
borders the hexagons seem conditioned the orientations the 
grains. Furthermore, possible detect the existence “‘macrosegrega- 
substructure, which forms enriched net solute closely linked with 
the boundary the grain. corresponds size the original grain and 
within there whole number Smialowski cells. 

The joint utilization the method etch pits and that epitaxial layers 
makes possible establish: (1) distribution the dislocation impurity 
associations starting from the interface and (2) their relation with segregated 
areas. 

While the interface the etch pits are concentrated along the borders 
the cell; towards the inner grain they concentrate macrosegregated areas, 
the vertices the cells, and alignments determining the substructure 
the macromosaic. This substructure closely linked with the vertices the 
Smialowski cells and lesser degree with their walls. The interpretation 
the distribution corrosion figures made accordance with the theory 
the formation dislocations through microsegregations 
cation. 

The effect crystalline orientation important all the phenomena 
pointed out. For that reason are present engaged experiments with 
single crystals known orientation. 
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NOTES 


THE 3820-A 1A’ TRANSITION PROPYNAL 


The first electronic absorption band system propynal sets 
about 4200 (origin ca. 4143 and extremely weak (ca. m-atm 
vapor). The second, rather stronger, system extends from about 
3000 with electronic origin From their intensities, and 
analogy with the spectra other aldehydes (especially formaldehyde), 
seems natural identify the two systems with transitions respectively 
the triplet and singlet states the first excited electronic configuration. 
formaldehyde, the effect such excitation make the nuclear equilibrium 
configuration non-planar (Brand 1956; Robinson 1956) and seemed 
interest see whether such behavior retained substituted formaldehyde 
which the substituent group would expected able conjugate 
(in the broadest sense) with the carbonyl. 

The 3820-A system was observed low dispersion Howe and Goldstein 
(1958), who identified the probable electronic origin and tentatively assigned 
frequencies some upper-state fundamentals. have re-examined the 
spectrum comprehensively, high dispersion, with 6-meter Ebert grating 
spectrograph (King 1958) resolving 0.1 the second order, and have 
extended the measurements include the and isotopic 
species. The molecule very nearly symmetric top, and the relatively simple 
rotational fine structure the bands, both and types, well resolved 
(see Fig. 1). This allows the spectrum analyzed rather greater detail 
than perhaps usual polyatomic molecules this size, and the complete 
analysis will take time, this note presents some the preliminary results. 

The spectral origin, 26163.1 the normal isotope, marked 
strong perpendicular-type band shown its rotational J-structure 
polarized along the inertial c-axis: thus the electronic transition 
notation appropriate the ground state which known planar from 
its microwave spectrum (Costain and Morton 1959). There only one strong, 
long upper-state progression, also type bands, frequency identified 
with the C=O stretching vibration. Shorter type progressions occur 
other frequencies, and the difficulties often encountered distinguishing 
between bands originating from the zero vibrational level the ground state 
and “hot’’ bands can here avoided the use A-type combination 
differences obtained from and subband branches. addition weak short 
progressions (of most two members) are observed frequencies which 
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spectrum being studied King and Moule (personal communication) 
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Fic. The zero-zero absorption band near 3821 photographed the 
second order Ebert grating spectrograph. Upper: appearance the band around the 
band center. Lower: five successive P-type subbands (AK —1) placed side side, showing 
details the J-structure. 


(2,2) 
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the first one-quantum members are type hybrid bands involving 
non-totally symmetric (out-of-plane) vibrations. 

The analysis the excited state vibration frequencies which have been 
identified the normal molecule summarized the table, which also lists 


TABLE 


Excited state state 


Class (in-plane) 


(C=C stretch) 1945.7 2106° 
(formyl C—H bend) 1119.5 1389° 
(C—C stretch) 943.7 
(C—C=O bend) 506.7 614 
bend) 189.2 205.5 
Class (out-of-plane) 
vio (formyl C—H deformation) 462.0 981.3 
(skeletal deformation) 345.9 260.6 


Frequencies taken from measurements of central intensity maxima given to the 
nearest 

6From the infrared; very weak or absent in the ultraviolet spectrum. 

eUncorrected for Fermi resonance with (vs+v9) at 1311.5 cm™!. 

dUncorrected for Fermi resonance with 2v10 at 943.6 cm™!. 


the analogous ground-state frequencies (values from the electronic spectrum 
agree with those from the infrared (Brand and Watson 1960; King and Moule 
1961)). The vibrations which appear the electronic spectrum are seen 
all those involving motion around the formyl carbon atom and include only 
one the acetylenic vibrations. The most striking changes effect which 
decreases about the same proportion formaldehyde, and the two 
out-of-plane vibrations and Any evidence non-planarity the 
excited state might therefore sought the latter. However, the pro- 
gressions these vibrations not exceed most two quanta, might 
otherwise expected grounds change shape were 
involved. There some doubt about the value due Fermi resonance, 
but its probable value 943 indicates some positive anharmonicity 
which is, however, still small compared with effects observed the 
spectrum. can observed three quanta sequences, and behaves 
normally. Thus, from vibrational evidence, excited propynal planar very 
nearly planar the equilibrium configuration; the height any possible 
inversion barrier does not exceed small fraction the zero-point energy 

The rotational A-structure strongly degraded longer wavelengths and 
the J-structure very slightly degraded shorter wavelengths. Analysis 
straightforward, but has revealed number interesting perturbations. 
least two types Coriolis coupling have been well characterized. The one 
involves the ground-state frequencies and which were inaccess- 
ible the infrared, and takes form very similar that identified there 
and (Brand and Watson 1960). close 0.63. The other found 
(a’, 2520 lies only few above another level, 
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symmetry whose rotational constant differs such that its K-type 
energy levels rise more rapidly than, and overtake, those Crossing 
occurs the region where the levels mix equally, and the subbands 
single electronic band appear double. Rotational perturbations probably 
intermediate between these cases appear other bands. 

Finally, propynal shows exceptionally clear-cut predissociation seen 
abrupt break the subband structure one perpendicular band. occurs 
bands are present normally, but and branches are 
wholly missing. This resembles similar case recently described Clement 
and Ramsay (1961) HNO. However, because the smaller rotational con- 
stants, the predissociation limit can narrowed down much more closely 
the present case, for the energy difference between levels and 
only Another difference between HNO and propynal lies 
the dependence the predissociation vibrational co-ordinates. HNO, 
the limit observed the same energy matter which vibrational state 
the molecule in. propynal the limit clearly observed (as line-broadening) 
only Bands involving transitions other vibrational states much 
higher energies retain their rotational fine structure full. The predissociation 
limit gives upper limit for the ground-state dissociation energy 
(3.545 ev, 81.74 kcal/mole). From the extremely 
sharp onset predissociation conclude that the repulsive, predissociating 
potential curve must cut the curve nearly horizontally, and that the upper 
limit given above must very close the actual value. 
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THE IONIZATION POTENTIAL 
HERZBERG 


few ago the ionization potential CH; was determined from 
Rydberg series the vacuum ultraviolet. The result was 9.843 (Herzberg 
and Shoosmith 1956). When the first absorption bands were obtained 
the flash photolysis diazomethane near 1415 (Herzberg and Shoosmith 
1959), was suspected that they might represent the first member 
Rydberg series similar the 1500-A bands However, because 
the shift shorter wavelengths, proved much more difficult 
obtain higher members the Rydberg series. means grating, coated 
for high reflectivity the 1200-A region, has now been possible extend 
the spectrum shorter wavelengths and observe higher members 
the Rydberg series which make possible establish the ionization potential 

discussed elsewhere (Herzberg 1961) the 1415-A group consists 
features shorter wavelengths which, all probability, represent two 
transitions. the new plates three further groups the same type 
observed which together with the 1415-A group form clear 
Rydberg series. Fig. which shows the spectrum from 1324 1219 
the three groups are readily recognized. addition, lines 1329 and 
1277 and bands 1271 and are present, and, faintly, some 
bands which are very prominent the spectrum the residue the 
flash photolysis. would expected, the spacing the individual subbands 
each group decreases for the higher series members confirming that these 
subbands correspond different electronic transitions with the same principal 
quantum number. 

The first strong band each group listed the second column Table 
the third column the wave numbers obtained from the formula 


are given. seen that this formula gives very good representation 
the observed bands confirming the Rydberg character this series. Therefore 
the limit the series 83851 10.396 represents ionization potential 
the free radical. 


TABLE 
Rydberg series 


3 70634¢ 70621 

76553 76562 

79241 79243 
80677 


80688 


“Band origin. 
*Issued N.R.C. No. 6502. 
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There remains the question whether not the ionization potential 
obtained the lowest ionization potential From the electron con- 
figuration the upper state the 1415-A band appears very likely that 
the main bands the Rydberg series correspond while the 
other features correspond and 2pm This interpretation 
agreement with the small size the Rydberg correction. this basis the 
state the ion which the series converges the *II state 
the ion which expected the ground state the ion linear. 
state linear molecule must expect Renner-Teller splitting. this 
splitting were large, for example (Dressler and Ramsay 1959), 
would mean that one component the state would have potential 
minimum for non-linear configuration, and this minimum would course 
below the minimum the linear component. The limit the Rydberg 
series must correspond the linear component since only one strong 
band occurs each member, that is, according the 
principle, the geometrical configuration the molecule the Rydberg states 
must the same the ground state. Therefore, there were strong 
interaction the true ionization potential would 
somewhat smaller than the value given. 

There are strong (but perhaps not entirely conclusive) reasons against the 
assumption that the ground state non-linear and that therefore the 
true ionization potential less than the value given. 

(a) Transitions from the linear ground state those components the 
Rydberg states which the molecule bent should give rise long pro- 
gressions the bending vibrations. such progressions were found; but 
this may have been due the low intensity these bands. 

(b) bent its ground state one would expect Rydberg series 
bent has been shown previously (Herzberg 1961) bent 
the lowest singlet state which can observed absorption very 
short time intervals after the beginning photolysis. Under these conditions 
Rydberg series singlet was found. Again this could due low 
intensity but would difficult see why singlet Rydberg series should 
much weaker than triplet Rydberg series. 

(c) Jordan and Longuet-Higgins (1961) the basis molecular orbital 
calculations similar those which have given excellent representation 
the low excited states have concluded that should linear. 
This also suggested the fact that the first excited singlet state 
and the first excited doublet state (a; the mole- 
cule linear spite the fact that according Walsh (1953) the 

tends bend the molecule while neutral this respect. 
other words single electron the bending orbital (as the ground state 
not sufficient make the molecule bent. 

the basis these considerations appears almost certain that the 
value given above represents the first ionization potential This value 
may compared with the electron collision value Langer, Hipple, and 
Stevenson (1954), who gave ev. The large discrepancy must 
ascribed the fact that the electron collision experiments formed 
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excited state with excitation 1.4 ev. There are several expected 
low-lying states from the configuration these the state 
probably the lowest and might well lie 1.4 ev. 
the new value for the ionization potential combined with the 
appearance potential from (viz. 15.39 ev) obtained inde- 
pendently Langer, Hipple, and Stevenson (1954) and Waldron (1954), 
one obtains immediately for the dissociation energy the value 
(or 113 kcal). However, again one must consider the possibility that 
formed excited state and therefore the value just 
given can only considered upper limit for the true value. From 
and the heat formation CH, together with and D(CH) one 
can immediately derive lower limit for namely kcal. 


greatly indebted Mr. Shoosmith for obtaining the plates 
which this discussion based, and Drs. Lossing and Douglas 
for useful criticisms. 
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LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 
previous the first day the month issue. proof will sent the authors. 


Magnetization Nutations Maser Oscillation Pulses 


The amplitude-modulated structure sometimes observed oscillation pulses from two-level 
solid-state masers has been investigated theoretically several reports. two these (Kemp 
1959; Yariv 1960), has been fairly well established that this modulation can attributed 
large-amplitude nutation the magnetization vector this the case, would 
reasonable expect the r-f. paramagnetic susceptibility remain pure imaginary 
qué antity (at resonance) for the duration the Thus, the case maser oscillation 
resonant cavity, the susceptibility component x’’ along the imaginary axis would either 
negative positive, depending whether the specimen was emitting power 
to, receiving power from, the cavity. the longitudinal component 
volume magnetization, equal the steady state. 

This situation self-evident the model proposed Kemp (1959). not evident, 
however, (1960) relationships. This paper concerned with examination the 
latter relationships for disclosure this feature. 

(1960) equations are, Gaussian units, 


(4) 


components the transverse magnetization and r-f. magnetic field, respectively. 
the electron gyromagnetic ratio, the sample filling factor, and the loaded 

the cavity. the cavity resonance and paramagnetic transitiona! frequency 
case. the spin-spin relaxation time for the paramagnetic specimen. The 

Mo, could evaluated from knowledge the F-center concentration, 
splitting-field and operational temperature. Note that perturbation term (such 
H,) necessary the system depart from state unstable equilibrium. Selected 


values were: 


5X10~ gauss, 
gauss, 
104, 
0.1, 


where the last term corresponds negligible interaction. 

The five time-dependent variables were evaluated the author with the aid Bendix 
G-15 digital computer. The (1951) method was employed. The variables 
were evaluated intervals seconds over range seconds, Values 
for and were then obtained (where M?, H,? together 
with the phase angle between and vectors. was found 90° phase shift was 
indeed preserved between the two variables and with leading for and 
lagging for The time-dependent behavior and was similar that 
shown Yariv, with the exception that only one zero-crossing the component was 
found for the parameters employed. addition, the peaks were displaced very slightly 
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from those (At seconds). The spin trajectories rotating co-ordinate frame 
would thus agreement with those discussed Kemp (1960) for weak spin-spin inter- 
action. fact, when equations were also programmed the computer, the derived 
time-dependent magnetization and r-f. magnetic field patterns were similar shape and 
duration those shown Yariv, for similar boundary conditions and circuit parameters. 


The author wishes thank Mr. Firth, Department Natural Philosophy, St. Andrews 
University, Scotland, for many valuable suggestions this topic. 


Git, X. 1951. Proc. Cambridge Phil. Soc. 47, 96. 
Kemp, J. C. 1959. J. Appl. Phys. 30, 1451. 
Yariv, A. 1960. J. Appl. Phys. 31, 740. 
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Observation Noctilucent Clouds Western Canada 


Only one authentic report noctilucent clouds over Canada known the writers 
(Vestine 1934). This may compared with the reported occurrences each year 
the European part the U.S.S.R. (Grishin 1957) and the systematic studies the luminous 
clouds made the Scandinavian Peninsula (Witt 1957). examination the 16-mm, all-sky 
camera records aurora over West-central Canada, taken during the failed disclose 
cloud structures this type. Since there much distortion the portions the photographs 
close the horizon, difficult distinguish between luminous cloud and faint aurora; 
and this itself not conclusive proof that the clouds were not present many occasions. 

Possible reasons for the lack observations noctilucent clouds over Canada were dis- 
cussed seminar early July. Some the reasons advanced were: 

(a) failure make observations times and directions most favorable their detection; 

(6) confusion between the luminous clouds and weak aurora; 

significant differences the upper atmospheric circulations between Northern Canada 
and Europe. 

The clouds are observed most frequently when the sun from 18° below the horizon, 
altitudes below 6°, and latitudes north 50° during the summer months. The low altitudes 
for maximum visibility, and directions toward the sun and the zone maximum auroral 
occurrence for Canada, make easy for inexperienced observer designate all faint 
luminous structures auroral origin. notable difference the regional occurrence the 
clouds, actually exists, might help explain their origin. 

Probably result these discussions, display the noctilucent clouds was observed 
the morning July 24. number photographs were taken. The clouds had the tenuous 
and shiny white appearance described European observers. The solar depression during 
the period observation varied between and 11°. The major portions the visible clouds 
were below angular altitude and azimuths toward the sun. 

This observation together with the special circumstances which induced 
search for noctilucent clouds suggest that they appear with much greater frequency over Northern 
Canada than the number reported occurrences suggest. While now too late the year 
for systematic observational programs started, suggested that such programs should 
initiated future years. 


GRISHIN. 1957. Izvest. Akad. Nauk S.S.S.R., Moscow. 
VESTINE. 1934. J. Roy. Astron. Soc. 28, 249. 
Witt. 1957. Tellus, 9, 341. 


RECEIVED AuGusT 9, 1961. K. D. BAKER 
INSTITUTE OF UPPER ATMOSPHERIC PHYSICS, B. W. CurRIE 
Puysics DEPARTMENT, 

UNIVERSITY OF SASKATCHEWAN, 

SASKATOON, SASKATCHEWAN. 


Can. J. Phys. Vol. 39 (1961) 


3 
j 
4 
| 
| 


NOTES CONTRIBUTORS 
Canadian Journal Physics 


MANUSCRIPTS 


General.—Manuscripts, English French, should typewritten, double spaced, 
paper The original and one copy are submitted. Tables and captions 
for the figures should placed the end the manuscript. Every sheet the manuscript 
should numbered. Style, arrangement, spelling, and abbreviations should conform the 
usage recent numbers this journal. Greek letters unusual signs should written plainly 
explained marginal notes. Characters set boldface type should indicated 
wavy line below each character. Superscripts and subscripts must legible and carefully 
placed. Manuscripts and illustrations should carefully checked before they are submitted. 
Authors will charged for unnecessary deviations from the usual format and for changes 
made the proof that are considered excessive unnecessary. 

abstract not more than about 200 words, indicating the scope the work 
and the principal findings, required, except Notes. 

References.—References should listed alphabetically authors’ names, un- 
numbered, and typed after the text. The form the citations should that used current 
issues this journal; references papers periodicals, titles should not given and 
only initial page numbers are required. The names periodicals should abbreviated 
the form given the most recent List Periodicals Abstracted Chemical Abstracts. All 
citations should checked with the original articles and each one referred the text 
the authors’ names and the year. 

Tables.—Tables should numbered roman numerals and each table referred the 
text. Titles should always given but should brief; column headings should brief and 
descriptive matter the tables confined minimum. Vertical rules should not used. 
Numerous small tables should avoided. 


ILLUSTRATIONS 


General.—All figures (including each figure the plates) should numbered con- 
secutively from up, arabic numerals, and each figure referred the text. The author’s 
name, title the paper, and figure number should written the lower left corner the 
sheets which the illustrations appear. Captions should not written the illustrations. 

Line drawings.—Drawings should carefully made with India ink white drawing 
paper, blue tracing linen, co-ordinate paper ruled blue only; any co-ordinate lines that 
are appear the reproduction should ruled black ink. Paper ruled green, yellow, 
red should not used. All lines must sufficient thickness reproduce well. Decimal 
points, periods, and stippled dots must solid black circles large enough reduced 
necessary. Letters and numerals should neatly made, preferably with stencil (do NOT 
use typewriting) and such size that the smallest lettering will not less than 
high when the figure reduced suitable size. Many drawings are made too large; originals 
should not more than times the size the desired reproduction. Whenever possible 
two more drawings should grouped reduce the number cuts required. such groups 
drawings, large drawings, full use the space available should made; the ratio 
height width should conform that journal page but allowance must 
made for the captions. The original drawings and one set clear copies (e.g. small 
photographs) are submitted. 

Photographs.—Prints should made glossy paper, with strong contrasts. They 
should trimmed that essential features only are shown and mounted carefully, with 
rubber cement, white cardboard, with space between those arranged groups. 
mounting, full use the space available should made. Photographs are sub- 
mitted duplicate; they are reproduced groups one set should mounted, 
the duplicate set unmounted. 


REPRINTS 


total 100 reprints each paper, without covers, are supplied free. Additional reprints, 
with without covers, may purchased the time publication. 

Charges for reprints are based the number printed pages, which may calculated 
approximately multiplying 0.6 the number manuscript pages (double-spaced type- 
written sheets, and including the space occupied illustrations. Prices and 
instructions for ordering reprints are sent out with the galley proof. 


Contents 


Fejer—The effects energetic trapped particles magnetospheric 
Magarvey and Roy Bishop—Transition ranges for three-dimensional 


Hector and Forsyth—On the statistical properties optical and radio 


Swanson and Friedberg—The electrical resistivity Mn;ZnC 
Axford and unifying theory high-latitude geophysical 
Cloutier, Bruce Flann, and Jane Manson—Study the effect 
Katzman—Changes the diurnal hour maximum the cosmic-ray 
Loh—Uncharged conducting toroidal ring uniform electric 
Heraldo Biloni—Substructures and dislocations produced during solidification 


Notes: 


Letters the Editor: 


Baker and Currie—An observation noctilucent clouds 


Printed Canada University Toronto Press 


1397 


1409 


1418 


1423 


1429 


1433 


1465 


1477 
1486 
1495 


1501 


1508 
1511 


1514 


1515 


aa 

% 


